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SURVEY OF HEAT TRANSFER TO NEAR-CRITICAL FLUIDS 
by Robert C. Hendricks, Robert J. Simoneau, and Ray V. Smith* 


Lewis Research Center 
SUMMARY 

This report surveys near-critical heat transfer, including the following areas: 
thermal properties; heat-transfer boundaries of the near-critical region; free, natural, 
and forced convection experiments; oscillations; geometric effects; parameters which 
appear to be significant to heat transfer in the critical region; and theories which have 
been proposed for the region. 

The report points out that the state of the art is such that completely satisfactory 
theories or correlations are not available for heat transfer in the critical region. It 
does, however, indicate the most promising and widely used correlations and suggests 
procedures for approaching a critical-point heat-transfer problem. 


INTRODUCTION 

Transport processes, particularly heat transfer, in the near-critical region have 
been of interest for about the last 15 years. Current emphasis stems from applications 
which require the use of a fluid in the near-critical condition, from inadequate informa- 
tion to produce satisfactory design expressions, and from an inadequate understanding 
of the mechanics which produce the peculiar behavior in the near-critical region. 

Applications which are current or proposed include the use of near-critical helium 
to cool the coils of superconducting electromagnets and superconducting electronic or 
power-transmission equipment, the use of supercritical hydrogen as a working fluid or 
fuel for both chemical and nuclear rockets; the use of supercritical water in electricity 
generating plants, and the use of methane as a coolant and fuel for the supersonic trans- 
port. 


*National Bureau of Standards, Cryogenics Division, Boulder, Colorado; work 
partially supported by U. S. Atomic Energy Commission. 



Before and during World War II, Ernst Schmidt and his associates conducted re- 
search in heat transfer by natural convection. They found the heat-transfer coefficient 
for a near-critical fluid to be quite high, because the specific heat and the compressi- 
bility both appeared to approach infinity at the critical point. The application they had in 
mind was the cooling of turbine blades in the then new jet engines. 

In the mid- 1950’s, a substantial number of papers began to appear which reported 
a rather wide range of results and proposed mechanisms underlying some of the unique 
heat-transfer behavior. Several authors suggested that the heat- transfer characteristics 
of a fluid near its critical point resembled boiling. This idea persists in the literature 
today, but the extent and influence of this boiling-like behavior still remains to be finally 
established. The dominant near-critical working fluids have been hydrogen, water, 
carbon dioxide, Freons, and ammonia. 

The reported research may be divided into two broad classes: forced and natural 
convection in heated tubes, and free or natural convection from heated wires and flat 
plates. The first class provided experimental turbulent-heat-transfer data for use in 
design and in establishing the reliability of correlations and theoretical analyses. The 
second class of studies focused more attention on the detailed mechanisms. The data 
from both these studies showed substantial contradiction, but there is an ever- increasing 
understanding of the conditions likely to produce such contradictions. 

Theoretical heat-transfer analyses have not been very successful; however, this 
approach still shows promise in unravelling the mysteries of the near- critical region. 

By contrast, flow oscillations, characteristic of this region, have been examined ana- 
lytically and found to be reasonably well predicted by using rather conventional concepts 
in mechanics. 

This report surveys near-critical heat transfer, including the following areas; 
thermal properties; heat-transfer boundaries of the near-critical region; free, natural, 
and forced convection experiments; oscillations; geometric effects; parameters which 
appear to be significant to heat transfer in the critical region; and theories which have 
been proposed for the region. 

The report is written with a twofold purpose: The first is to bring the information 
in this region together and to examine these data from a designer’s point of view. The 
second is to examine the flow mechanisms underlying the behavior. Within the limita- 
tions of knowledge in this field, this examination lends itself not so much to selection of 
prime variables influencing the flow, but, more to recommendations for future work in 
the area. 

Since this is a survey report, the authors are indebted to the entire community of 
researchers in the fields of near-critical heat transfer and near-critical thermody- 
namics for the material presented herein. In particular we are indebted to several in- 
dividuals who made available their original figures for use in this report. They are 
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Dr. J. M. H. Levelt Sengers, National Bureau of Standards; Dr. J. V. Sengers, the 
University of Maryland; Dr. R. S. Brokaw, NASA Lewis Research Center; Dr. R. S. 
Thurston of Los Alamos Scientific Laboratory; Dr. P. Giffith, Massachusetts Institute 
of Technology; Dr. M. Cumo, Centro Studi Nucleari, Rome, Italy; Dr. C. Gazley, 

Rand Corp; Dr. G. E. Tanger, Auburn University; Dr. J. R. Bartlit, Los Alamos 
Scientific Laboratory; Mr. S. S. Papell, NASA Lewis Research Center; Dr. L. W. B. 
Hall and J. D. Jackson, University of Manchester, Manchester, England; Mr. R. G. 
Deissler, NASA Lewis Research Center; Dr. N. Zuber, New York University; Dr. A. J. 
Cornelius and J. D. Parker, Oklahoma State University; Dr. R. C. Williamson, NASA 
Electronics Research Center. 

NEAR-CRITICAL FLUID PROPERTIES 
Thermodynamics of the Critical Point 

The exceptional heat-transfer behavior of a near-critical fluid must ultimately be 
due to the influence of the unusual property behavior of a fluid near its critical point, as 
shown in figure 1. This is manifested both in the changes of the thermal properties 
themselves, and in the modifications of the flow structure caused by these changes. For 
any analysis to succeed, a good knowledge of the thermal properties is required. 

The classical approach to the thermodynamics of the critical point is well docu- 
mented in Hirschfelder, Curtiss, and Bird (ref. 1) and in Rowlinson (ref. 2). For ex- 
cellent surveys of some of the more recent ideas, the reader is referred to Hammel 
(ref. 3), Widom (ref. 4), and Sengers and Levelt Sengers (ref. 5). 

The Van der Waals equation of state 


P = - a lP 2 (1) 

1 - bjp 

where a^ is an attraction force constant and b^ is a repulsion force constant, is 
representative of the classical approach. (All symbols are defined in appendix A. ) 
Among other things, at the critical point the Van der Waals equation (eq. (1)) yields an in- 
finite specific heat at constant pressure and an infinite thermal expansion coefficient 
(-3 In p/3 T)p, which Schmidt (ref. 6) cited as reasons why the critical fluid was an at- 
tractive heat-transfer medium. 

The Van der Waals equation of state (eq. (1)) yields a cubic curve for the critical 
isotherm on the pressure-density (P-p) plane and a quadratic for the coexistence curve. 
In reality (ref. 3), the critical isotherm is more nearly fifth degree, and the coexistence 
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Figure 1. - Typical thermodynamic and transport properties of near-critical parahydrogen. Ratio ot pressure to critical pressure, 1.05. 
(Data taken from refs. 10 and 17). 


curve is almost a cubic. This does not affect the primary condition for determining the 
critical point, that of equal saturation densities. It does, however, interject uncertainty 
as to the nature of the higher derivatives of P with respect to p. 

The Van der Waals equation also yields a finite discontinuity for the specific heat at 
constant volume C y at the critical point. Experimental evidence by Bagatskii, et al. 
(ref. 7) indicates that C y has a logarithmic infinity with a superimposed finite discon- 
tinuity along the critical isochore. Widom’s recent qualitative survey of his activities 
(ref. 4) proposes a three-dimensional lattice gas model to account for these difficulties. 

Since near-critical thermodynamic data are somewhat scarce, the heat-transfer in- 
vestigator is likely to have to extrapolate from other regions or other fluids. This 
prompts a comment on the law of corresponding states. To a good first-order approxi- 
mation, it can be said that the thermal properties of all fluids collapse to a single curve 
when normalized by the appropriate critical constants. There are some problems. The 
second-order deviations are more severe near the critical point. For spherical non- 
polar gases such as hydrogen and nitrogen, Z c = 0. 304 and 0. 292, respectively; where- 
as, for the polar molecule water, Z Q = 0. 224, approximately a 25 percent difference. 
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The hydrocarbons fall in between (ref. 1). When dealing in corresponding states, it is 
best to subdivide the fluids by molecular type. For quantum fluids, like helium, corre- 
sponding states should be avoided. 

A final area of importance to heat transfer is the nonequilibrium phenomenon asso- 
ciated with the critical point. Maass (ref. 8) found that a hysteresis loop in density near 
the critical point could be obtained by first heating and then cooling along an isobar, as 
shown typically in figure 2. The results were very stable and reproducible even under 


Density 


Temperature 

Figure 2. -Typical near-critical isobar exhibiting hysteresis loop 
according to Maass 1 experiment (ref. 8). 

conditions of stirring. Another hysteresis example can be found in light scattering. 

Light scattering caused by severe density fluctuations near the critical point causes the 
fluid to become opaque as it passes through the critical point. The growth of the opaque 
condition, known as critical opalescence, is quite different depending on whether heating 
or cooling through the critical point is occurring. (For an excellent visual record of the 
critical opalescence phenomenon, the reader is referred to a motion picture produced 
by Siemens Aktiengesellschaft, Postfach, West Germany, entitled "Boiling and Evapora- 
tion Phenomena with Water. ”) These experiments demonstrate that, near the critical 
point, large relaxation times are required for a thermally disturbed system to return to 
equilibrium. One consequence of this is that the near-critical fluid under the dynamic 
conditions of heat transfer can be expected to experience some degree of thermodynamic 
nonequilibrium. This, in turn, will result in some uncertainty in applying an equation 
of state. In a private communication, Walter G. Zinman (Grumman Aircraft) argues 
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that fluctuations are inherent to the thermodynamic critical point and that therefore the 
experimentalist may not be able to precisely locate such a point. 

Also, since the paths are different, the heat-transfer process involved in cooling a 
near- critical fluid may be quite different than that associated with heating the same fluid. 

Most of the truly severe behavior, such as the singularities in C p , C y , X, and /3j, 
occur precisely on the critical point. Normally, operation is not precisely on the criti- 
cal point. This will tend to attenuate the influence of some of the anomalies. On the 
other hand, the overall large property changes, the suggested first-order effect, persist 
over a considerable region near the critical point. 


Pressure-Density-Temperature Data - Equations of State 

The availability of actual P-p-T data near the critical point varies considerably 
from fluid to fluid. Probably the most detailed investigations near the critical point have 
been in carbon dioxide by Michels, Blaisse, and Michels (ref. 9) and in hydrogen by 
Goodwin, et al. (ref. 10). For other fluids of interest, like nitrogen, there is little, if 
any, actual data. For the most part, the heat-transfer researcher must rely on prop- 
erties computed from an equation of state. 

Equations of state come in all sizes and shapes. The most common are the virial 

type 


P = A(T)p + B(T)p 2 + C(T)p 3 + D(T)p 4 +. . . 


( 2 ) 


Obert (ref. 11) and Hirschfelder, et al. (ref. 1) have good surveys. They list coeffi- 
cients for various equations and fluids; however, they were not determined with the crit- 
ical point specifically in mind. In cryogenic fluids, the Benedict-Webb-Rubin (ref. 12) 
equation of state, as modified by Strobridge (ref. 13) and Roder and Goodwin (ref. 14), 
has been very popular: 


n 3 n 4 n 5\ 2 


P = (fiTp -+ [ tfln jT + n 2 + + — + — p* + 

T ^,2 ,p4/ 


n 9 n 10 n ll\ ~ n 16 p ' 


(RngT + Urj + j + 

rp2 m 3 rpl 


+ ngTp + 


i2 n 13 n 14\ " n 16 p 

— + + |e 

T 2 t 3 rp4 


,5 6 ,n\ 

p + n 15 p (3) 
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TABLE I. - SOURCES OF NEAR-CRITICAL PROPERTIES 


Properties recommended as computable by source authors 

Thermal 

conductivity, 

X 

X 

X 

X 

X 




X 

X 



X 

Viscosity, 

V 

X 

X 

X 





X 

X 



X 

Entropy, 

S 

X 

X 




X 


X 


X 

X 


Expansion 

coeffi- 

cient, 

P \3Ty p 


X 


1 


X 






X 

Specific 
heat at 
constant 
volume, 

C v 


X 


X 


X 


X 


X 



Specific 
heat at 
constant 
pressure, 

c p 


X 


X 


X 


X 





Enthalpy, 

H 

X 

X 

X 



X 


X 


X 

X 


Pressure- 

density- 

temperature, 

P-p-T 

X 

X 

X 

X 

X 

X 


X 


X 

X 

X 

Source reference 

. .. . 

Transport 

properties 




33,35, 

176 

e- 

CM 

1 

1 

1 

• 

I 

1 

i 

i 

i 

i 

i 

i 

i 

33, 35 

24,25, 

30,33, 

35 

t 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

28,33, 
35, 185 

Equation 

of 

state 

o 

c— 

tH 

eg 

eT 

eg 
~ c- 

CO r-t 
„ 1— » 

cT 

rH tH 

CO 

c- 

co 

lO 

r- 

tF 

m 

t- 

c— 

oo 

-Q 

179 to 182 

CM 

CM 

nj 

oT 

183 

184 

185 

Critical constants 

Density, 

Pc’ 

g/cm 3 

eg 

CD 

o 

© 

co 

o 

o 

tH 

CO 

o’ 

CO 

m 

CO 

o* 

0. 1625 

0.210 

i 

i 

t 

i 

« 

i 

0.322 

0.468 

0. 557 

0. 582 

0.236 

Temperature, 

T c- 

K 


1 

126.2 

154.8 

190.7 

370.0 

i 

i 

i 

i 

i 

i 

647. 3 

304.2 

385. 1 

418.7 

406. 1 

Pressure, 

P c- 

atm 

1 

9 

33. 5 



50.14 

45.66 

42. 1 

i 

i 

i 

i 

i 

i 

CM 

CO 

CM 

72.90 

40.69 

31.45 

112. 3 

Fluid 

1 

S 

3 

’q; 

X 

c 

<D 

ho 

O 

U 

T3 

>. 

X 

Nitrogen 

Oxygen 

Methane 

Propane 

Kerosene 

Water 

Carbon 

dioxide 

w 

CM 

1 

C 

o 

<L> 

U 

{*< 

Freon- 114 

Ammonia 


7 


Limited range of application. 

’No known equation of state. 

Decomposition potentialities near the critical point must be recognized when using this fluid. 














An alternate form of the Benedict-Webb-Rubin (BWR) equation was used by Vennix 
(ref. 15) to correlate his near-critical-methane data. The coefficients are of a different 
type than equation (3), and the exponentials are linear in p and T. The main difficulty 
with these or any curve- fit equation is determining other properties, such as specific 
heat C , since this requires derivatives which are often not very satisfactory near the 
critical point. For specific discussions of accuracy, the reader is referred to the in- 
dividual references (table I). 

Table I lists source references for equations of state for the fluids discussed here- 
in. The table includes an indication as to whether the particular authors suggest com- 
putation of other useful thermodynamic properties. In some cases, P-p-T are not used 
as primary variables. The table also lists transport properties. These will be taken 
up in the next section. For some fluids, such as hydrogen, entire computerized prop- 
erty packages are available (e.g. , refs. 16 and 17 and a National Bureau of Standards 
(NBS) property package based on refs. 10 and 14). 

Recently there has been considerable interest in the use of the Ising model for 
ferromagnets to describe the lattice gas near the critical point (refs. 18 and 19). This 
is discussed in the surveys (refs. 3 to 5) mentioned earlier. In drawing the analogy 
with the Ising model, chemical potential appears to be a more fundamental variable than 
pressure. Green and his co-workers (refs. 20 and 21) have offered the following equa- 
tion, in which chemical potential is expressed as a function of density and temperature: 

Jp(p,T) - ju(P c > t )]p c 

i i "^2^ 1 

= l x ll lh ( x l) ( 4 ) 

1 



/r 

where t = (T - T c )/T c and Xj = tf [(p - P c )/p c J 

In reference 21 the function h(xj) had yet to be determined analytically but could be 
seen graphically (fig. 3(a)). In a more recent paper (refs. 22 and 23), Vicentini- 
Missoni, Levelt Sengers, and Green have succeeded in developing a closed-form ex- 
pression for h(Xj) 
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H T = 29.929 DEG C 



(a) Graphical representation of h(x,l (b) Analytic correlation of h(x,) for cartx>n dioxide (ref. 22). 

(ref. 21). 


Figure 3. - Scale-reduced chemical potential as function of scale-reduced density. (Data reproduced from refs. 21 and 22). 


The results are shown for carbon dioxide in figure 3(b). Thus a closed-form equation of 
state is now available which is directed specifically at the critical region and which is 
subject to a very few adjustable constants, /3g, 6^, Ej, Eg, and x q . Of these, only x q 
is very sensitive. Its major limitation is range. It is presently limited to ±30 percent 
in A p and -0. 01 < t < 0. 03. Thus in a practical property package, such as references 
16 and 17, equations (4) and (4a) will have to be merged with the more common virial 
equations such as equation (3). 


Transport Properties 

Most of the remarks made concerning P-p-T data can be made about transport 
property data. Again, the more popular fluids are better documented. There are 
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near-critical thermal conductivity data for carbon dioxide by Sengers (ref. 24) and by 
Guildner (ref. 25), for hydrogen by Differ and Roder (ref. 26), for methane and ammonia 
by Sokolova and Golubev (refs. 27 and 28), and for nitrogen by Ziebland and Burton 
(ref. 29). Viscosity has been measured for carbon dioxide by Kestin, Whitelaw, and 
Zien (ref. 30) and for hydrogen by Diller (ref. 31). 

Prior to Sengers' work, it was thought the excess of viscosity or thermal conduc- 
tivity above its "atmospheric" value at a given temperature and pressure was a function 
of density alone: 


X(P,T) - X*(p - 0,T) = fj(p) 

(5) 

?](P,T) - T?*(p - 0,T) = f 2 (p) 

(6) 


where fj(p) and f g(p) may be determined to almost the critical temperature for a gas 
by a density variation of the form: 


fj(p) — + ^2 ^ P + ^2^ + • • • (5a) 

f 2 (p) = ’ll P + v' 2 ln P + V 2 P 2 + • • • (6a) 

Of course, this required a precise knowledge of p as a function of P and T; never- 
theless, it allowed a very simple representation for viscosity and thermal conductivity. 
Furthermore, it could be reduced to general form by the law of corresponding states. 
This is still true away from the critical point, except possibly for quantum liquids such 
as hydrogen and helium (see refs. 26 and 32). The functional relations of equations 
(5) and (6) have been established for most fluids of interest by Thodos and co-workers 
and are summarized in a paper by Stiel and Thodos (ref. 33). Equations (5) and (6) re- 
main valuable as a base line. In a private communication, H. J. Hanley of NBS, 
Boulder, Colorado, examined the hydrogen viscosity data of Diller (ref. 31) and inter- 
preted the results to mean that there exists, at best, a weak anomaly; however, it does 
not occur at the critical density. Therefore, for viscosity the best evidence indicates 
only a weak anomaly, which can be disregarded for most fluid flow and heat-transfer 
calculations. 

All the thermal conductivity data cited above exibit an "anomalous spike" in the 
conductivity at the critical density along near-critical isotherms, as shown in figure 4 
(taken from ref. 34). 
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Figure 4. - Thermal conductivity of carbon dioxide exhibiting the 
anomalous spike. (Data reproduced from ref. 34). 


The computation of the anomaly is still in developmental stages. Brokaw (ref. 35) 
suggests treating a near-critical fluid as a dissociating polymer and that conductivity 
consists of two parts 


X - + X r (7) 

Here X r represents the contribution due to diffusion of the dissociating clusters and 
X f the normal conductivity expressed by equation (5). The results of Brokaw's analysis 
are 

« 

A r^ D ^( C p) r < 8 > 

where D is the self-diffusion coefficient, D^ n is the binary diffusion coefficient of the 
hypothetical polymer, and is the specific heat in excess of the low-pressure value. 
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Density, p, g/cm 3 

Figure 5. - Calculated and experimental thermal con- 
ductivities of carbon dioxide. (From ref. 35.) 


Brokaw's calculations are compared to Sengers' data in figure 5. Unfortunately, the 
theory requires a good equation of state for computing and Dj n /D. So far it has 

not been tried for other fluids. 

Even though the data of Diller and Roder (ref. 26) for hydrogen and the data of Hanley 
and Childs (ref. 32) for helium place in considerable doubt the dependence of A - X q 
on p alone, even away from the critical point, equation (5) will probably have to be 
treated as the best available in all cases, except for liquid hydrogen and helium. 

For the present, we can recommend computation of thermal conductivity by equa- 
tion (7) using equation (5) to obtain X f and by Brokaw's theory, where an adequate equa- 
tion of state is available, for X . Viscosity can be computed directly from equation (6). 

Many additional related properties (vapor pressure, heat of vaporization, surface 
tension, and others) can be found in the NBS compendium by Johnson (ref. 36). 
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Pseudocritical Properties 


The unusual property behavior in the near- critical region shown on figure 1 can be 
examined as an extension of saturation properties. 

Frequent reference is made in near-critical work to the transposed critical tem- 
perature T*. This is normally defined as the temperature where the specific heat C p 
attains a maximum for a given supercritical pressure (see fig. 1). The locus of the 
T* points for hydrogen is shown in figure 6 and seems to be a very natural extension of 
the saturation curve. It is not, however, the same as an extension of the vapor pressure 
curve by the Clapeyron equation using enthalpy difference for latent heat (fig. 6) 

(refs. 37 and 38). Near the critical point, however, the T* locus can be considered a 
pseudosaturation temperature. The T* locus can be considered to approximately sep- 
arate the pseudoliquid from the pseudovapor. Notice also on figure 6 that the T* (or 
max C ) effect has pretty well diminished at a P/P c corresponding to the maximum in 
the Joule-Thorn pson inversion locus. Thermodynamically then, it can be said that the 
critical region is confined within the inversion locus. 

In order to use pseudoboiling models, certain two-phase quantities in addition to the 
vapor pressure curve must be defined. Typical of these is latent heat of vaporization. 



Figure 6. - Thermodynamic region of influence of the critical point on heat transfer. Basic locus is for hydrogen. (Clapeyron ex- 
tension of saturation curve from ref. 37.) 
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30 40 50 
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Figure 7, - Pseudo-two-phase properties; extrapolation procedure tor equivalent latent 
enthalpy of vaporization. (From ref. 40.) 


Figure 7 shows a way of defining a pseudoheat of vaporization from the enthalpy curves, 
as used by Thruston and others (refs. 39 and 40). Other pseudo-two- phase quantities, 
such as ’’saturated” liquid and vapor densities, can be approximated in the same way 
(refs. 41 and 42). 

In summary we can say: 

(1) Classical Van der Waals models do not accurately describe the state of a near- 
critical fluid, but may be satisfactory in first-order considerations. 

(2) The discrepancies in the law of corresponding states are accentuated in the 
critical region, and the law should be used only between fluids of similar molecular 
structure. 

(3) A thermally disturbed near-critical fluid requires long relaxation times to re- 
turn to equilibrium. 

(4) Equations of state of the conventional virial type can be made to accurately map 
a near-critical P-p-T surface; however, they are not good for thermal properties 
requiring derivatives of the equation of state. 

(5) A new analytic equation of state involving chemical potential, density, and tem- 
perature, based on the Ising lattice gas model, has been developed for the near-critical 
region. 
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(6) The evidence for the existence of a strong anomalous spike in thermal conduc- 
tivity at the critical point is overwhelming. The spike can be estimated by treating the 
fluid as a dissociating polymer and superimposing the results on the conventional 

X - X* = i^(p) curve (eq. ( 5 )). 

( 7 ) If a similar spike occurs in viscosity, it is sufficiently weak to be ignored. 
Thus, 77 - 77* = fgOo) (eq. (6)) will adequately describe viscosity. 

(8) Pseudo-two- phase thermal properties can be generated by extrapolation of 
saturation properties and relations into the near-critical region. 


HEAT-TRANSFER REGIONS 

The heat-transfer problem will be broken into four parts: gases (regions I and la), 
liquids (region II), and two phase (region IE), which are the three regions adjacent to 
the critical region; and region IV, the near-critical region itself. These four regimes 
are illustrated in figure 8. Although this map is based on extensive hydrogen data, a 



Figure 8. - Heat-transfer regions as function of inlet con- 
ditions. Region I, gas or fluid; region la, gas (low 
temperature); region II, liquid (fluid); region III, two 
phase,- region IV, near critical. 
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similar map could be anticipated for other fluids. In order to determine the heat trans- 
fer and pressure drop, the region of interest must first be determined. Such a procedure 
is discussed and flow-charted in appendix B. In the following sections, before going to 
the near-critical heat-transfer problem (region IV), recommendations are offered for 
the adjacent regions I, la, n, and m. 


Region I - Gas-Fluid 


The data of McCarthy and Wolf (refs. 43 and 44) and Taylor (ref. 45) for hydrogen 
and helium give extensive heat-transfer and pressure-drop coverage of region I. The 
recommended equation is that of Taylor, which includes entrance effects and reduces to 
the successful equation of McCarthy and Wolf (refs. 43 and 44) at X/d > 25. 


/T x-(0. 57 - 1.59 d/X) 
Nu = (Nu) b (-^\ 

T b/ 


( 9 ) 


where 


(Nu) b = 0.023(Re)°' 8 (Pr)°- 4 (10) 

Taylor (ref. 46) and McCarthy and Wolf (ref. 44) also determined friction factor for 
the data of region I. The equation of Taylor is recommended. It approaches the classic 
equation of Koo-Drew and McAdams at low T w /T b - Taylor’s equation is 


i_ 

2 


0.0007 + 


0.0625 


(Re) 


0. 32 
w 



(ID 


Region la - Near-Saturation, Low -Temperature Gas 

This region represents an extension of region I, with gas-like behavior and strong 
property variations. Perroud and Rebiere (ref. 47) and J. Gladden of Lewis (personal 
communication) investigated pressures ranging from atmospheric to several atmos- 
pheres, with inlet temperatures near saturation. Perroud’s equation is 


16 


(Nu) b = 0.0032(Re)°- 95 (Pr)°- 4 

Taylor was able to correlate much of these data by using equation (9). 


Region II - Liquids 

Region II of the map (fig. 8) would be classified as liquid, with a process at low- 
to- moderate temperature differences. The conventional techniques of handling liquid 
heat transfer as given by McAdams (ref. 48), and also by equations (9) and (10), are 
expected to be valid (perhaps modifying the constant 0. 023). 


Region III - Two Phase 

The two-phase regime, region HI in figure 8, has been actively explored for many 
cryogens and noncryogens. Perroud and Rebiere (ref. 47), Hsu, et al. (ref. 49), 
Giarratano and Smith (ref. 50), Hendricks, et al. (refs. 41 and 51), and others, have 
explored cryogens. Noncryogen works are extensively covered in a text by Tong 
(ref. 52). 

Heat transfer to two-phase cryogenic fluids at high heat fluxes can be adequately 
correlated by the theoretical approach of Hsu, et al. (ref. 49) and the semiempirical 
approach of Perroud and Rebiere (ref. 47). The latter is quiet similar to that of 
Hendricks, et al. (refs. 41 and 51). In most cases, the pressure drop is due to momen- 
tum changes, and friction pressure drop is difficult to assess. 

Hsu, suggested that two-phase heat transfer could be treated with traditional single- 
phase correlations if the thermal and transport properties were weighted according to 
the void fraction; that is, 


(Nu)^= 0.023(Re)“- 8 (Pr)“- 4 


( 12 ) 


where 


cp = ct l(Pl + a x <p x 


and cp represents any and all the fluid properties in equation (12). The reference con- 
dition x was determined empirically for hydrogen and, so far, has not been checked 
with other fluids (see ref. 49). 
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The equations given by Perroud have been applied to tubes and channels, as well as 
to several different cryogens. In all cases, the basic form remains the same: 

( Nu W Sr ) 1/3 = f (X tt ) = a + — L - (13) 

b + cxt t 

however, the function f(x tt ) is altered slightly by differences in fluids and geometry 
(e. g. , for hydrogen flowing through a tube, a = 1 . 16, b = 0. 16, c = 0. 3, and dj = 1; and 
for hydrogen flowing through a channel, a = 1. 08, b = 0. 135, c = 0. 20, and dj = 1. 0 
(ref. 53). 

Phenomenologically, at sufficiently high T w /T b values, vapor can be expected at 
the wall and fluid particles (filaments) in the core, while just the opposite is true for 
low T w /T b values. The occurrence of large temperature differences or large T w / T b 
is more natural in cryogenic fluids, and is one of the major reasons for the once pre- 
valent opinion that heat transfer to cryogens and noncryogens differed. In both boiling 
and near-critical fluids, the range in T w - T b or T w /T b is pertinent to heat transfer 
and will be considered more fully in the section NEAR-CRITICAL HEAT-TRANSFER 
REGION. 


Boundaries of Region IV - Near-Critical Region 

It is difficult to define the boundaries which separate the near-critical region 
(region IV) from its adjacent regions. There are several reasons for this: (1) For most 
fluids, the data are not sufficient; (2) the transition is not abrupt and sharp demarcations 
are difficult to determine; and (3) the near-critical influence will persist further into the 
adjacent regions, depending on the path (process) the fluid took to arrive at a given state 
point (see appendix A). This last point is particularly elusive. A gas can be precooled 
at critical pressure down to, and possibly even below, the transposed critical tempera- 
ture T*; and, when run in a heated-tube experiment, it will behave as a precooled gas 
(i. e. , similar to region I). On the other hand, if liquid is heated up through region IV 
to the same state point in a heated tube, the results will be quite different. Part of the 
near-critical heat-transfer problem is how the fluid arrives at a given state point (i. e. , 
prior history does make a difference). This is discussed more fully in appendix A. 

Nevertheless, classification of the heat-transfer regions by state conditions is con- 
venient, uesful, and generally reasonable. For hydrogen (as shown in fig. 8), the pres- 
sure boundaries are 0. 8 < P/P c < 3. The lowest temperature boundary is T ga ^, cor- 
responding to P/P c = 0. 8, and the upper boundary is the vicinity of T*, the transposed 
critical temperature. Other fluids have not been explored extensively enough to establish 
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their boundaries or to confirm the universality of the hydrogen boundaries. In the ab- 
sence of data, the hydrogen boundaries can probably be taken as reasonable for other 
fluids. 

A comment is in order here concerning boundaries. What is meant in establishing 
these boundaries is that outside of them the influence of the critical point can be con- 
sidered to be negligible and conventional variable property correlations will prevail. It 
does not mean that in every case within region IV conventional methods will fail. Some 
combinations of parameters within region IV will be amenable to conventional approaches; 
but, in general, this region will require analyses directed specifically at the near- 
critical heat-transfer phenomenon. 

The problem of heat-transfer regions can be summarized in the following manner: 

(1) The near-critical heat-transfer region and its adjacent regions can be delineated 
in an approximate manner on a temperature- entropy state diagram. 

(2) It is not possible to set precise near-critical heat-transfer boundaries for the 
following reasons: 

(a) Data are generally insufficient. 

(b) The transition from region to region is not abrupt. 

(c) Most important, the extent of influence of the critical point on heat transfer 
is a function of the process (or path) by which the fluid approaches the critical point. 

(3) Heat transfer in the adjacent regions can be correlated by conventional methods 
appropriate to that region. In fact, the correlatable nature of the adjacent regions helps 
determine the near- critical boundaries. 

NEAR-CRITICAL HEAT-TRANSFER REGION 
Pecularities of the Near-Critical Region 

In region IV (fig. 8), standard techniques of correlating data break down. The ordi- 
dary Dittus-Boelter equation (eq. (10)) does not correlate near-critical forced convection 
phenomena; ordinary Rayleigh relations do not correlate pool data; standard boiling 
equations exhibit discontinuities; and oscillations are commonplace. Standard theoreti- 
cal techniques are no better than the empirical techniques. The problem, at the risk of 
oversimplification, is that the heat-transfer coefficient has a strong and complex tem- 
perature dependence unlike an ordinary gas (fig. 9). 

The early experiments of Schmidt, Eckert, and Grigull (ref. 6) found free and 
natural convection in the vicinity of the critical point to exhibit a sharp increase in heat- 
transfer coefficient. On the other hand, Powell (ref. 54) reported a sharp minimum in 
heat-transfer coefficient in the T* region for forced convection flow of liquid oxygen 
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(a) Forced convection. 



(b) Free convection - poo! boiling. 


Figure 9. - Comparison between subcritical and supercritical heat -transfer behavior in forced and free convection 
(refs. 41, 51, 59, 62, and 631. 


and nitrogen. Since then, many other investigators have found similar "peaks” in the 
axial wall temperature profile (refs. 42, 55 to 57). In direct opposition, there are 
several researchers (refs. 58 to 61) who have reported a maximum heat-transfer coef- 
ficient (fig. 9(a)). The forced convection results of figure 9(a) are for hydrogen (ref. 42) 
and carbon-dioxide (ref. 59), while figure 9(b) represents pool heat transfer for hydro- 
gen (ref. 62) and water (ref. 63). The similarities inwall temperature profiles for sub- 
and supercritical hydrogen data are also shown. 

There are two associated results which may shed some light on this problem and 
which bear further investigation. First, both Shitsman (ref. 56) and Yamagata, et al. 
(ref. 57) showed heat-transfer coefficients, when plotted against bulk enthalpy, to be a 
minimum near the critical enthalpy and also reported pressure oscillations in the same 
region. Second, all the reports of a maximum in heat-transfer coefficient came from 
experiments (refs. 58, 60, 61, and 64 to 66) in which the temperature difference be- 
tween the wall and the bulk fluid was small when compared to experiments in which a 
minimum occurred (refs. 56 , 57 , 67, and 68). 

Hsu (ref. 69) suggested that the two results can be thought of qualitatively in boiling 
terms. When the temperature difference is small, it can be likened to nucleate boiling, 
a region of very good heat transfer, thus the maximum. When the temperature differ- 
ence is large, it can be compared to film boiling, a region of poor heat transfer, thus 
the minimum. In support of this line of thinking, examine Hauptmann’s data (ref. 70) 
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Heat-transfer coefficient, h, W/(crrn(K) 



Figure 10. - Effect of free-stream temperature and wall heat flux on near- 
critical heat-transfer coefficient for carbon dioxide. Ratio of pressure 
to critical pressure, 1.025; free-stream velocity, 0.45 meter per second. 
{Data taken from ref. 70.) 



.S2 150 \— 


100 \— 



1500 2000 2500 

Bulk enthalpy, kJ/kg 

Figure 11. - Variation of heat -transfer coefficient with heat flux in 
critical region exhibiting both a maximum and a minimum. Mass 
flow rate, 700 kilograms per square meter per second. (From 
ref. 74.) 
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shown in figure 10. The small-temperature-difference data (i. e. , usually small heat 
flux q) exhibit a clear maximum, while the large-temperature-difference data (q = 8. 4 
W/cm 2 ) show, if anything, a minimum. A similar trend was found by Styrikovich, et al. 
(refs. 71 and 72) and is shown in figure 11. Wood (ref. 66) showed the heat-transfer 
coefficient to be a maximum near the critical point when plotted as a function of bulk 
temperature, and Kahn (ref. 73) showed the maximum to be a function of wall tempera- 
ture. Shiralkar and Griffith (ref. 74) have found the heat-transfer coefficient to be 
strongly dependent on heat flux when wall temperatures are above T* and bulk tem- 
peratures below T*. 

Pressure oscillations are a natural phenomenon of this regime and at times can be 
quite large, 0. 3 P test or 400 psi (27. 2 atm) in nitrogen tetroxide as reported by 
McCarthy, et al. (ref. 75). Hines and Wolf (ref. 76) reported pressure oscillations of 
sufficient amplitude to damage their heated tubes. 

Pressure drops are large and are primarily due to momentum changes (ref. 41). 
Friction losses become increasingly important as the heat-transfer process moves away 
from the near- critical regime. 

The similarities in heat transfer, pressure drop, wall temperature profiles, and 
pressure oscillations between subcritical and the near-critical regimes is remarkable. 

In the early phases of the hydrogen work, the sounds of two-phase flow and supercritical 
flow were recorded (ref. 77). The sounds were similar, but supercritical sounds were 
not as "noisy. " Goldman (ref. 78) has noted the same phenomenon. Researchers in 
this area (refs. 71, 79, and 80) often think of the fluid as pseudo-two- phase. Some of 
the strongest evidence in support of a pseudo-two- phase fluid comes from figure 6, from 
visual studies (refs. 62, and 81 to 83), and from heat-flux-against-temperature plots 
such as figure 9 (refs. 62, 63, and 83). Griffith and Sabersky (ref. 81) and Knapp and 
Sabersky (ref. 82) found the heated globules to be easily taken for bubbles. The film 
supplement to the work of Graham, et al. (ref. 62) gives the viewer a statistical feel 
for the nature of these similarities, as opposed to the instantaneous picture provided by 
the published photographs (e. g. , ref. 82). Cumo, et al. (ref. 84) have documented the 
the change in fluid structure for Freon from atmospheric to the near-critical region as 
shown in figure 12. Nishikawa and Miyabe (ref. 83) published a set of photographs for 
the nucleate-film-boiling cycle and a comparative set at the reduced pressure of 
P/P c = 1. 065. The similarities at the film-boiling heat fluxes are particularly striking. 
The heat-flux-against-temperature-difference plots cited above, particularly those of 
Holt and Grosh (ref. 63), show the heat-transfer coefficient to be a function of heat 
flux - a phenomenon not found in normal gas heat transfer, but a trademark of boiling. 
One possible explanation of the unusual wall temperature profiles discussed earlier is a 
"boiling" model. 

The strongest argument against such a pseudo-two-phase fluid is thermodynamic; 
a supercritical fluid in equilibrium is clearly single phase. Pseudo-two-phase-fluid 





Reduced pressure, P/P c : 

Temperature, t, c C: 

Mass flow per unit area, G, g/cnV 
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Figure 12. - Heat transfer to Freon-114 at various 




advocates recognize this, and most, such as Miyabe (ref. 83) and others, claim that it 
is a nonequilibrium situation. Experimentally it is very interesting that the single-phase 
proponents use the same visual evidence presented by the two-phase advocates. 
Hauptmann (ref. 70) uses his own very dramatic photographs and reexamines the other 
visual experiments (refs. 70 and 81), and concludes that all of the unusual results can 
be explained in single-phase terms. A careful examination of all the available photo- 
graphs (refs. 62, 70, and 81 to 83) reveals as many nonsimilarities to boiling as there 
are similarities. Really, the analyst has his choice because a strong case can be made 
for either model; however, it is clear that there are no sharp discontinuities in the 
heat-transfer data(refs. 67 and 63). Any analytic model (pseudo- two- phase or single- 
phase) should work throughout region IV (fig. 8) and provide smooth transition from 
liquid to gas if properly formulated. 

The heat-transfer results change from apparatus to apparatus. For example, fig- 
ure 13 illustrates changes in heat transfer for hydrogen as tube diameter is changed. 

At the larger diameter, the L/d effects are large and not accounted for by the pseudo- 
fluid or single-phase models. Thus the results appear to be system-dependent, so 
geometry effects must also be carefully assessed. 



(a) Hess and Kunz correlation dependence on pseudo- 
Martinelli parameter and diameter {ref. 41). 



(b) Taylor equation behavior for two different tube 
diameters. 


Figure 13. - Some geometry effects in near-critical - 
fluid heat transfer. 
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In summary we can say 

(1) Some characteristics of the near-critical region are 

(a) Nonlinearities in heat flux against temperature difference 

(b) Wall temperature excursions (spikes) 

(c) Similarities to the two- phase regime 

(d) Oscillations 

(e) Large momentum pressure drops 

(f) System-dependent results 

(g) Failure of conventional correlations 

(2) Modeling should be required to work throughout region IV, and should exhibit a 
smooth transition to the adjacent regions shown in figure 8. 

Heat Transfer in Free Convection Systems 

The free convection pool studies discussed in this section deal primarily with heat 
transfer from small test sections, such as wires, filaments, and plates. In a subsequent 
section, the heat transfer from natural convection loops will be considered. 

The effect of the wide variations of properties appears to have a more direct in- 
fluence in the pool-heat-transfer case. This is probably true because in this situation 
temperature differences are usually small. Pool results almost universally show en- 
hancement in heat transfer near the critical point. Some typical experimental data where 
enhancement occurred are those of Dubrovina and Skripov (ref. 85) shown in figure 14. 


20xl0 3 



Figure 14. - Free convection heat transfer for carbon dioxide. Data show heat -transfer-coefficient behavior as the properties 
near the critical point, when temperature differences are very small. (From ref. 85). 


25 



TABLE II. - FREE 


Investigator 

Date 

of 

work 

Refer- 

ence 

Fluid 

Geometry 

Doughty and Drake 

1956 

186 

Freon- 12 
(CC1 2 f 2 ) 

Horizontal wire 

Boniila and Sigel 

1960 

187 

n -pentane 
< C 5 H 12> 

Horizontal plate 

Griffith and 
Sabersky 



Freon-114A 

(C 2 Cl2F 4 ) 

Horizontal wire 

. 

Skripov and 
Potashev 

1962 

188 


Vertical wire 

Brodowicz and 
Bialokoz 

1963 

189 

Freon- 12 

(cci 2 f 2 ) 

Vertical plate 

Holt and Grosh 

1963 

■ 

h 2 o 

Horizontal wire 
and vertical 
ribbon 

Fritsch and Grosh 



h 2 o 

Vertical ribbon 

Dubrovina and 
Skripov 

1964 

85 

co 2 

Horizontal wire, 
horizontal and 
vertical ribbon 

Knapp and 
Sabersky 

1965 



Horizontal and 
vertical wire 

Nishikawa and 
Miyabe 

1965 

83 

O 

Horizontal wire 

Graham, Hendricks, 
and Ehlers 

1965 

62 

H 2 

Horizontal 

ribbon 

Larson and 
Schoenhals 

1966 

192 

h 2 ° 

Vertical ribbon 

Goldstein and Aung 

1967 

86 

co 2 

Horizontal wire 

Daniels and 
Bramall 

1965 

193 

co 2 

Horizontal 
Monel wire 
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- POOLS 


Range of conditions 

Constants 

Ratio of pressure 
to critical 
pressure, 

p/p c 

Ratio of temper- 
ature to critical 
temperature, 

T/Tc 



0.28 to 1.7 

0.77 to 1.04 

n - 1/3 
C = 0. 14 

0.72 to 1.5 

0.65 to 0.74 

C = 0. 

172 


(fluid) 

n - 1/3 


0.71 to 1.06 
(surface) 

01 AT= ln(p b /p w ) 

0.42 to 1. 15 

0.725 to 1.08 



0.82 to 1.25 

1.00 to 1.02 



1 

0.98 to 1.07 



1.5X10" 4 to 1.25 

0.75 to 1.0 
(fluid) 

0. 75 to 2.0 
(surface) 



1.03 to 1.06 




0.99 to 1.3 

1.00 to 1.03 



1.03 to 1.40 

0.93 to 1.09 



1.00 to 1.33 

0.98 to 1.63 

Pseudocritical 
properties analysis 

0.32 to 1.38 

0.76 to 1.19 

0.25 <n 

£0.35 

1.01 to 1.03 

0.98 to 1.04 

a = 1/3 ^ 
b - 0.247 
c = 0. 137 
x - 0. 5 ^ 

► Eq. (15) 

1.03 to 1.21 

0.94 to 1.09 




0. 79 to 1. 16 


0.86 to 1.03 



























The enhancement of the heat-transfer coefficient is very pronounced but also limited to 
a very small temperature range. Also the enhancement peaks somewhat away from the 
critical point and in the region of the transposed critical temperature. 

Unlike a forced convection system, the free convection system has no constraining 
boundaries. Consequently, it is able to respond to favorable property variations with 
enhanced heat transfer. 

Although boiling-like behavior has been reported for pools (refs. 63 and 81), the on- 
set of this behavior was not characterized by a marked change in the heat-transfer coef- 
ficient, as is the case with the onset of nucleate boiling. This can be seen in the super- 
critical curves of figure 9. Although the mechanisms for nucleate boiling and near- 
critical heat transfer may be similar, some aspects of the energy transport must be dif- 
ferent. 

Some interesting near-critical free convection experiments are those of references 
82, 86, and 87 using a small-diameter horizontal wire. Knapp and Sabersky (ref. 82) 
were the first to photograph an oscillation back and forth between laminar and turbulent 
flow while heating near-critical carbon dioxide. Their data indicated a sharp increase 
in heat transfer in going from the all- laminar to the all-turbulent region; however, the 
transition through the oscillating region appeared to be smooth. The authors indicate 
that this is probably because the data were average data, and that they suspect an oscil- 
lation in wire temperature to go with the laminar-turbulent oscillation. Goldstein and 
Aung (ref. 86) reported similar oscillations but no sharp increase in heat transfer. The 
phenomenon remains unexplained at present; however, such oscillations may well be 
fundamental to our understanding of the oscillations noted in forced convection and natural 
convection loops. 

Recently, Hahne (ref. 88) conducted the classic free convection experiment between 
two horizontal plates, using supercritical carbon dioxide. He found "spindle- like" 
structures, rather than the classic Benard cells. 

The pool experiments in the near-critical region are summarized in table II. Most 
of the data analyses of the papers shown in table n employ rather conventional dimen- 
sionless groups. In some cases, there have been modifications, primarily to account 
for the variable properties in the boundary region (refs. 89 to 91). There appears to be 
general agreement that somewhat away from the critical and transposed critical points 
the conventional correlations will hold for all fluids and for all systems, as discussed 
in the section HEAT-TRANSFER REGIONS. Thus, away from the critical point, 

McAdams basic equation (ref. 48) 


(Nu) f = C(Ra) n 


(14) 


should work. Closer to the critical point, modifications of the following basic form are 
employed: 
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Specific details are listed in table II. Very close to the critical point, Simon and Eckert 
(ref. 92) recommend 



In summary, 

(1) Pool results almost universally show an enhancement in heat transfer near the 
critical point. In all cases, the temperature differences are small and the enhancement 
occurs over a limited temperature range and peaks somewhat away from the critical 
point. Large temperature differences diminish the peak. 

(2) Free convection systems have no constraining boundaries and consequently re- 
spond favorably to property variations. 

(3) Boiling-like behavior in near-critical pools is not accompained by the marked 
increase in heat-transfer coefficient of nucleate boiling. 

(4) Most researchers employed conventional techniques away from the critical point. 
However, close to the critical point, modifications such as equations (15) and (16) are 
employed (see table II). 

(5) Laminar-to-turbulent oscillations have been observed, but the origin of this 
phenomenon has not been explained. These oscillations may well determine the character 
of near-critical oscillations found in natural and forced convection systems. 


Heat Transfer in Loops - Natural Convection Systems 

Historically, there has been considerable interest in the use of a natural convection 
loop or column (fig. 15). Here a fluid in the supercritical region can operate at high 
heat flux. 

Experiments and analyses in natural convection loops begin with the work of Schmidt 
(refs. 6, and 93 to 95). As in pools, the results of loop experiments are universally 
an enhancement of heat transfer near the critical point when compared to noncritical 
fluids under similar heat-transfer conditions. And as in forced convection systems, 
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oscillations appear to be dependent on the proximity to the critical point, as will be dis- 
cussed in a later section. The laminar-to-turbulent-flow oscillations noted in near- 
critical pools could be the origin of the oscillations noted in natural convection loops. 

Further evidence of similarities to forced convection are the body force effects at 
low flow (discussed under geometry effects) and heat transfer at low AT = T w - T^. 

In an effort to correlate near-critical data, Holman and Boggs (ref. 96) rearranged 
Schmidt’s basic loop equations to the form 


Nu = 16(Re) 2 (Pr)(Gr) 



(laminar) 


(17) 


Nu = 0. 079(Re) 11 / 4 (Pr)(Gr)' 1 |— V-\ (turbulent) 


(18) 


To test the geometry effects, Tanger, et al. (ref. 97) compared sulfur hexafluoride 
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Figure 16. - Correlation of near-critical sulfur hexafluoride data in natural convection 
loops of different geometries. (From ref. 97.) 


data obtained in two different loops. For a fixed and d, the geometry factors of 
equations (17) and (18) of loop 1 were approximately 1. 8 those of loop 2. They found no 
effect due to geometry varations. For the critical region, they recommend 


Nu(Gr)j 

Pr 


0. 00982(Re) 11//4 


(19) 


The results are shown in figure 16. Merlini (ref. 98) found a very strong influence of 
the transposed critical temperature on both heat transfer and oscillations. He used the 
experimentally determined film-to- center line fluid temperature ratio T x /T^ to corre- 
late his data 


where 


/T V 5 - 27 

Nu = Const(Re) 0, 65 (Pr) 0, 15 L* \ 


0.287 


Const = < 


T x 

— > 1.2 

T* 


0. 196 


T x 

— < 1.2 

T t 


(19a) 
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TABLE m. - NATURAL CONVECTION - LOOPS 


Investigator 

Date 

of 

work 

Refer- 

ence 

Fluid 

Geometry 

Range of conditions 


Ratio of 
pressure 
to critical 
pressure, 

p/p c 

Ratio of 
temperature 
to critical 
temperature, 
T/T c 

Schmidt, Eckert, 
and Grigull 

1939 

6 

Ammonia 

(NH 3 ) 

Loop: 0. 59 in. i.d. , 
7.4 ft long 

0.52 to 1.29 

0.6 to 1.25 

Schmidt 


93 

h 2 ° 

Turbine blades 
(hollow core) 

1 

0.77 to 1.07 

Van Putte and Grosh 


194 

h 2 o 

Test section: 5 ft long, 
0. 245 in. i.d. 

0.03 to 1.25 

Up to 1. 12 

Merlini 


98 

h 2 o 

Loop: 0.742 in. i.d., 
42 ft long 

1.06 to 1.24 

0.735 to 1.016 

Holman and Boggs 

1966 

96 


Loop: 0.43 in. i.d. 
49 ft long 

0.84 to 1.5 

0.88 to 1.39 

Harden and Boggs 

1964 

155 

Freon- 114 
(C 2 C1 2 F 4 ) 

Loop: 0.93 in. i. d. , 
49 ft long 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Cornelius and 
Parker 

1965 

151 

Freon- 114 
(C 2 C1 2 F 4 ) 

Loop: 0.93 in. i.d. , 
49 ft long 

1.0 to 1.1 

0.86 to 1.0 

Tanger, Lytle, and 
Vachon 

1968 

97 

Sulfur 

hexafluoride 

(SF 6 ) 

2 loops: 43 in. i. d. , 
28 and 38 ft long, 
varied vertical 
distribution between 
condition and heater 

0.98 to 1.64 

0.995 to 1.14 

Hahne 

1965 

195 

co 2 

Closed inclined tube 

Constant volume and temper- 
ature experiment 
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Because of this equation’s extreme dependence on T /T, , the reader should become 
familiar with reference 98 and exercise caution where applying equation (19a). The ex- 
periments in natural convection loops are summarized in table in. 

The results obtained in free convection pools bear many similarities to the results 
for natural convection loops; however, there are also many similarities between loops 
and forced convection systems. 

In summary, 

(1) In free convection pools and natural convection loops, close proximity of the fluid 
to the critical point appears to enhance heat transfer. This is possibly the result of a 
free boundary system being able to respond to enhancing property variations when the 
temperature differences are low. As will be shown, forced convection systems at low 
temperature differences exhibit similar trends. 

(2) Reasonable success has been achieved in correlating the data of both the pools 
and natural convection loops by variable property modifications to the standard equations 
for the given configuration. 

(3) The pool experiments have produced considerable support for a pseudo- two- phase 
modeling of a near-critical fluid. On the other hand, the success of single-phase corre- 
lations indicates either approach may be successful. 

(4) To date no single correlation has been tested against all the available data. 

(5) Since these experiments operate very close to the critical point, there is a con- 
tinuning need for improved property data and equations of state. 

(6) System oscillations which are commonplace in a natural convection loop near the 
critical point are similar to those observed in forced convection systems. 


Heat Transfer in Forced Convection Systems 

The experimental work in near-critical forced convection can be broken into two 
major categories. First are the conventional heated-tube experiments which have been 
used so successfully in determining Nusselt correlations for gases. In the second cate- 
gory are experiments which examine one or two details of the heat-transfer process in 
order to explain the mechanism of near-critical heat transfer. 

Heated-tube experiments . - In the past 17 years, over 30 experiments in forced con- 
vection heat transfer to supercritical fluids, which fall into the broad classification 
known as heated-tube experiments, have been conducted. All these experiments were 
very similar in design, operating procedure, measurements, and to some extent, re- 
sults. These experiments can be surveyed as a group, pointing out specifics where 
necessary. The range of these experiments, which by now is quite extensive, is sum- 
marized in table IV. The most popular fluids have been water, hydrogen, and carbon 
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TABLE IV. - FORCED CONVECTION - HEATED- TUBE EXPERIMENTS 


Investigator 

Date 

of 

work 

Refer- 

ence 

Fluid 

Range of conditions 

Correlation attempts 

Ratio of 
pressure 
to critical 
pressure, 

P/P c 

Ratio of bulk 
to critical 
temperature, 

V T c 

Ratio of wall 
to critical 
temperature, 
T /T 

V J c 

Heat flux, 

q ’ 2 
W/cni 

Mass flow 
per unit 
area, 

G, 

g/cm^-sec 

Character- 
istic length 
(diameter) 
d, 
cm 

Length- 

diameter 

ratio, 

L/d 

Equa- 

tion 

Constants 

Monroe, Bristow, and 
Newell 

1952 

80 

°2 

1.03 to 2.74 

0.99 to 1.29 

1.07 to 1.60 

0.29 to 12. 8 

2. 72 to 71.5 

12.7 

75 to 220 

None 



C ha If ant 
Randall 
Goldman 

1954 

1956 

1961 

196 

197 

198 

|h 2 o 

1.56 

0.74 to 1.25 

0.74 to 1.60 

32 to 940 

210 to 550 

1.27 to 1.90 

10 to 120 

None 



Dean and Thompson 

1956 

79 

N 2 

0.61 to 1.22 

0. 64 to 0. 66 

0. 8 to 9. 0 

65 to 810 


4.75 

Average 

None 




Bringer and Smith 

1957 

100 

co 2 

1.12 

0.97 to 1.06 

1.02 to 1. 12 

3. 2 to 32 

54 to 220 

4.57 

0 to 133 

(15) 

C = 0.0266, a = 0.77, 

■(*$«“• ■>-». 

X . ![<T* - T b )/(T W - T b )] 

Powell 

1957 

54 

°2 

1.20 to 1.44 

i 

0.63 to 3. 6 

2.5 to 7.4 

190 to 800 

300 to 2740 

4.93 

31 to 371 

Mixing 

length 



Miropoiskii, Picus, and 
Shits man 

Miropolskii and Shitsman 

1966 

1957 

105 

102 

j H 2° 

0.018 to 1.29 

0.43 to 1.07 


29 to 2320 

17 to 3000 

2.9 to 8.2 

10 to 40 

(17) 




Dickinson and Welch 

1958 

64 

H 2° 

1.09, 1.41 

0.58 to 1.25 

0.65 to 1.33 

870 to 1820 

220 to 340 

7.62 

25 to 180 

(15) 

C = 0.023, a = 0.8, b = 0.4, 
p = 0, x = 0 

Morra and Barger 
DelBene and Barger 

1958 

1959 

199 

200 

}cci 2 f 2 

0.98 to 1.30 

0.97 to 1.22 

0.98 to 1.72 

0.55 to 10.7 

21 to 71 

8.90 

0 to 131 

None 



Schmidt 

1951 

67 

h 2 ° 

0. 78 to 1.37 

0. 73 to 1.25 

0. 77 to 1.55 

29 to 82 

70, 120 

5.0, 8.0 


None 




Koppel and Smith 

1961 ' 

59 

| C °2 

1.000 to 1.03 

0.96 to 1.06 

1.10 to 2. 10 

6.3 to 63 

26 to 135 

4.93 

2.5 to 93 

None 



Petukhov, Krasnoschekov, 
and Protopopov 

1961 

60 

co ? 

1.20 to 1.46 

0.98 to 1.15 

1.01 to 1.16 

3.4 to 12.0 

68 

6.70 

5 to 97 

(19) 

' a - -0. 11, b = 0.35, c = 0, 
d 0. 35 

(Nu = Nu(Re, Pr, f» 

Thompson and Geery 

1961 

68 

p-H 2 

3.58 to 7. 10 

0.92 to 3.93 

2.0 to 26.6 

18 to 1300 


4.93 

20.6, 41.2 

(15) 

See ref. for constant 


Szetela 

, 1962 

143 

P -H 2 






6.58 to 11.4 

1 to 34 

None 



j 





Shits man 

1963 

56 

h 2 o 

1.05 to 1.15 

0.85 to 1.01 

| 0.91 to 1.35 

<115 

30 to 150 

8.0 

10 to 186 

(17) 

j 


Domin 

1963 

58 

h 2 o 

1.01 to 1.19 

1 si. 12 

<1.53 

1.15 to 115 

<516 

2.0, 4.0 

Average 

(17) 



Wood 

1963 

; 66 

co 2 

l 1.004, 1.03 

0.94 to 1.01 

0.99 to 1.08 

1 

3.2 to 20.5 

144 to 226 

22.9 

30.8 

Mixing 

length 



Krasnoshchekov and 
Protopopov 

Krasnoschekov, Protopov, 1 
Feng, and Kuraeva 

1964 

1966 

106 

201 

|co 2 

1.06, 1.33 

0.96 to 1.26 

0.97 to 2.6 

4.6 to 260 


4.08 

1 to 51 

(19) 

a = 0, b = 0, c = 0.3, 
d = f(T b , T w , T*) 

(Nu = Nu(Re, Pr, f)) 

Finn 

1964 

108 

C 3 H 8 

1.13 

0.84 to 1.07 

0.96 to 1.43 

3.15 to 38 

53 to 173 

7.65 

3. 5 to 67 

(19) 

i a = -0.02, b - 0.20, 
c -- -0.27, d = 0.54 

Miller, Seader, and 

Trebes j 

1965 

I 

55 

p-h 2 

2.13 to 13.3 

0.85 to 3.09 

,16 

223 to 3910 

250 to 1180 1 

5.34 

2.5 to 55 

(21) 

C = 0.0204, x - 0.4, 
a = 0. 00983 

Holman, Rea, and 
Howard 

1965 

202 

cci 2 f 2 

0.90 to 1.27 

0.95 to 1. Osj 

1.00 to 1.301 

1 J 

0.47 to 27 

17 to 44 

12. 7*44. 4 
annu lus 


(20) 

See ref. for constants 


Hendricks, Simoneau, 
and F riedman 

1 

1965 

203 

P-»2 

5.30 to 13.3 

1.00 to 3. 36 j 2.5 to 20.0 

58.5 to 1630 

1 

28 to 840 

5.35 to 11.1 

3.4 to 78.2 

(15) 

C = 0.021, a = 0. 8, 
b = 0. 4, x = 0.5, p = 0 

Swenson, Carver, and 1 

Kakarala 

1965 

61 

H 2° 

1.03 to 1.88 

0.54 to 1.31 1 0.57 to 1.42 

1 

20 to 182 

54 to 218 

9.43 



20.2 to 182 

: 

(19) 

1 

Slightly different form 
(see ref.) 
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TABLE IV. - Concluded. FORCED CONVECTION - HEATED- TUBE EXPERIMENTS 


Investigator 

Date 

of 

work 

Refer- 

ence 

Fluid 



Range of conditions 

I Correlation attempts 

Ratio of 
pressure 
to critical 
pressure. 

P P c 

Ratio of bulk 
to critical 
temperature. 
T b T r 

Ratio of wall 
to critical 
temperature. 
T w T c 

Heat flux, 

q. 

W cm 2 

Mass flow 
per unit 
area, 

G. 

g/cm 2 -sec 

Character- 
istic length 
(diameter) 
d, 
cm 

Length- 

diameter 

ratio. 

L d 

Equa- 

tion 

Constants 

Bishop, Krambeck, and 
Sandberg 

1965 

204 

h 2 ° 

0.75 to 0.975 

0.97 to 1.24 

1,02 to 1.41 

23 to 256 

66 to 330 

2.54, 5.08 

30 to 385 

(15) 

C 0.0073. a 0.886. 
b 0.61.x 0.5, p -- C . 

plus L/d term 

Hendricks, Graham, 
Hsu, and Friedman 

1966 

41 

p-h 2 

0.43 to 4.37 

0.75 to 8.38 

2.0 to 15.0 

33 to 540 

14 to 510 

4,75 to 12.8 

5 to 114 

(22) 


Tanaka, Nishiwaki, 
and Hi rata 

1967 

65 

co 2 

1.06 

0.987 to 1.02 
.997 to 1.03 


a l . 43 to 5.20 
b 1.43 to 8. 14 

50 to 100 

10.0 

5 lo 95 

Mixing 

length 


Yamagata, Nishikawa, 
Hasegawa, and Fujii 

1967 

57 

h 2 ° 

1. 15 

0.607 to 1.25 

0.915 to 1. 50 

11.6 to 93 

32.3 to 108 

10.0 

10 to 190 

(19) 

Constant 

f[(T> - T b )/(T W - T b )] 

(see ref. ) 

Schnurr 

1969 

101 

C° 2 

1.003 to 1.03 

0.97 to 1.04 

1.07 lo 1.39 

25 to 189 


2. 64 

1 to 75 

(15) 

C 0.0266. a 0.77. 

ws-wt". P-* 

X = IpT- - T b )/<T„ - T b )] 
(Nu • Nu(Re. Pr. f)> 

Armand, Tarasova, 
and Konkov 

1959 

205 

h 2 ° 

1.04 to 1. 19 



18. 7 to 73.5 

45 to 65 

6.0. 8.0 

3. 1 to 37 

(17) 

C 0.0215 

Shitsman 

1962 

206 

H 2° 

1.03 to 1. 19 

0.89 to 1.03 

0. 82 to 1.01 

b 29 to 580 

10 to 250 

8.0 

4 to 18 

(19) 

a -0.25. b 0.33. 
c 0. d 0.35 
(Nu Nu(Ito, Pr. f)) 

Krasnoschekov and 
Protopopov 

' 1959 

145 

c° 2 

1.23 to 1.37 





6.70 

5 lo 95 

(19) 

For regions see ref. 

Glickstein and Whitesides 

1966 

207 

ch 4 

0.3 lo 1.2 

0.3 to 1.8 


to 374 

89 to 362 

8.5 

99 

(22) 

Not >j t dependent 

Heech c 

1967 

208 

Keros t' nn 

0.83 to 3.0 

0.43 to 0.46 

0.44 to 1.17 

to 250 

250 to 265 

6.0 

51.5 


St(Pr) 0,6 0.0034 

Tube in a rectangular 
passage 

Kafengauz 

1969 

153 

Di isopropyl 
cyclohexane 

1.5 to 2.5 





0.8 to 2.3 

13 to 38 

None 







Von Berg, Williamson, 
and Edeskuty 

1969 

209 

n 2 

0.81 to 1.22 

0.70 to 1.23 


1.9 to 7.5 

1.1 to 6.8 

34. 8 

131 

(10) 


Hay 

1969 

210 

He 

1.7 to 10.2 

0.80 to 1.50 

0.85 to 1.60 


3.9 to 22 

0.25 

600 

None 


Puvarnin and Malkina 

1969 

211 

Various 

hydro- 

carbons 

2 to 4 



35 

400 to 3200 

1 to 3 

100 

(20) 


Touba and McFadden 

1966 

212 

H 2° 

1 . 006 to 1 . 06 

0.916 to 1.01 

0.926 to 1.11 

55.5 to 165 

115 to 525 

0.622 and 
0.94 

212 and 
151 


f" -1 — 1(„°- 2 0.0068 
K-vJ 

-expp.ta^ 1 - - o.aoijj 

a Heating. 

^Cooling. 

Also in an unpublished English Ministry of Technology Report by Beech and Ziebland entitled "Heat Transfer to K< 

jrosine at Sup 

ercritical and 

Subcritical Pressures. " 
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dioxide; and the flow is normally vertically upward. The emphasis in many of these ex- 
periments was on data correlation. 

The correlation attempts can be treated in a general manner as to form. Early work 
centered on reference temperature modifications to the basic Dittus-Boelter equation: 

(Nu) x = C(Re)*(Pr)k^ (20) 

where x means that the fluid properties are evaluated at T x = T^ + x(T w - T^). Eckert 
(ref. 99) suggested that x may be a function of a dimensionless temperature frequently 
called the Eckert parameter 



( 21 ) 


Bringer and Smith (ref. 100) and Schnurr (ref. 101) have plotted these functions, as 
shown in figure 17. While this treatment handles quite a bit of data, it has never been 
generalized to other data and no analysis is available to predict the function. 



Figure 17. - Effect of relative proximity to the critical point on Eckert reference tem- 
perature (refs. 59, 101, and 135). 
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Another early correlation directed specifically at heat transfer near the critical 
point was that of Miropolskii and Shitsman (ref. 102); using water experiments they pro- 
posed 


(Nu) b = 0.023(Re)“- 8 (Pr)^t 


( 22 ) 


where (Pr) . is the minimum Prandtl number, defined as 
mm 


(Pr'min H 


(Pr) b (Pr) b < (Pr) w 

(Pr) w (Pr) w < (Pr) b 


(22a) 


This amazingly simple modification to the Dittus-Boelter equation produced some re- 
markable results. Miropolskii and Shitsman (ref. 103) made a plot (redrawn in fig. 18) 
for water at 230 atmospheres (P c = 218. 2 atm) using equation (22). Figure 18 shows a 
sharp peak in heat-transfer coefficient at — T* and the peak decreases with in- 
creased heat flux until, with an order of magnitude increase in q, there is practically 
no peak. Compare this to the data from Hauptmann's experiment (ref. 70) for carbon 
dioxide flowing over a flat plate, shown in figure 10. The qualitative similarities be- 
tween these data and the Miropolskii and Shitsman calculations (fig. 18) are striking. In 



Figure 18. - Calculations of heat-transfer coefficient for 
water by Mi ropolsky and Shitsman correlation (eq. (22)). 
Ratio of pressure to critical pressure, 1.052; mass flow 
rate, 100 grams per square centimeter per second; tube 
diameter, 30 millimeters. (Data from ref. 103.) 
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Hauptmann's data, however, the heat-transfer coefficient first increased, then decreased 
with heat flux. The data also showed a marked transverse body-force effect which, of 
course, the correlation could not comprehend. Shitsman (ref. 104) successfully ex- 
tended this correlation to include the water data of Dickinson and Welch (ref. 64), the 
oxygen data of Powell (ref. 54), and the carbon dioxide data of Bringer and Smith 
(ref. 100). Particularly interesting was the fact that the correlation comprehended both 
the maximum in heat-transfer coefficient reported by Dickinson and Welch and the mini- 
mum reported by Powell. Later, however, using water, Shitsman (refs. 103 and 105) 
also discovered experimentally the same temperature "spikes, " shown typically in fig- 
ure 11, that Powell (ref. 54) reported for oxygen. He called this a region of "impair- 
ment" to heat transfer and likened it to boiling crisis. He cautioned his readers that 
equation (22) would not correlate this region despite the fact that it had earlier handled 
Powell's similar conditions. The authors have tried this equation on the hydrogen data 
of Hendricks, et al. (ref. 41) with unsuccessful results. 

Another correlation directed at the critical point is of the form 



where (Nu)^ is the Dittus-Boelter value (eq. (10)) and Cp = (H w - Hj 5 )/(T w - T^). Ver- 
sions of equation (23) have been used for carbon dioxide (refs. 60 and 106), for water 
(refs. 61 and 107) and for propane (ref. 108). The key characteristic would seem to be 
the integrated specific heat Cp. 

The concept of using not only integrated specific heat, but all physical properties on 
an integrated average basis was put forth by Brokaw (ref. 109) for a reacting sys- 

tem. The properties would be expressed as 



where cp( t) is any fluid property. This has a tendency to smooth out the sharp near- 
critical property changes, as shown in figure 19 for density. Since near-critical fluids 
can be thought of as dissociating, this approach could have merit. Furthermore, if the 
temperature profile is known, or a reasonable assumption for the profile is made, the 
integral average property becomes 


(24a) 


^i = - fo <?[ T (y)] dy 

r 
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Temperature, T, K 

Figure 19. - Integral average properties as defined by Brokaw (ref. 109) 
in comparison to equilibrium properties, for parahydrogen. 


Hess and Kunz (ref. 110) proposed a correlation of the form 


(Nu) x 


= C(Re)°' 8 (Pr)°' 4 ' 


1 + a 



(25) 


This has worked reasonably well for hydrogen (ref. 55). 

Most authors have been reluctant to propose "pseudoboiling” models for forced con- 
vection heat transfer. Hendricks, et al. (refs. 41 and 42) have suggested modification 
of two boiling concepts, the Martinelli parameter (ref. Ill) and the Stermann parameter 
(sometimes called boiling number) (ref. 112), for use in the supercritical region: 


(Nu) f =(Nu) fm f( Xtt ’,Sr) 


(26) 




where 


and 



(26a) 


Sr=-SL^- 

a b AH 


(26b) 


where the subscript 1 refers to the liquid or heavy species. While the functional re- 
lation has not been specified, these parameters have shown some success in grouping the 
data. 

Wilson (ref. 113) used a somewhat similar concept: 


1 

(Nu) w 


= a + bN 


sk 


(26c) 


where 


and 


N _ Q - ^ - H inlet> 
sk d>(H g -H z ) 

a = 1. 033xl0“ 4 ; b = -1.708x10" 3 N gk < - 0. 2 
a = 7. 868X10" 4 ; b = 1. 613xl0" 3 N gk > - 0. 2 


(26d) 


If an equation is listed as a correlating equation in table IV, it is intended to imply 
that the particular author used a nondimensional correlation which is one of the basic 
forms of equations (20) to (26), as indicated in the table. The reader is referred to the 
actual reference for the specific form of the correlation. 

In forced convection, in addition to subdividing regions by state conditions as in 
figure 8, some success has been had by Yamagata, et al. (ref. 114) and others (refs. 41 
and 42) in determining correlating regions by flow parameters. Figure 20 shows a typi- 
cal map of heat flux against mass flow. Basically, for combinations of high flow or low 
heat flux, conventional correlations of the type shown in equation (9) will work even in 
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Figure 20. - Typical map of heat flux as a function of mass flux. 
(From refs. 74 and 114.) 


region IV of figure 8. For high heat flux or low flow, however, conventional correlations 
fail and the peculiarities of the near-critical region, such as wall temperature spikes 
appear. This information helps considerably in reducing the size of the region of non- 
correlation (region IV). Unfortunately, for the present, only water and carbon dioxide 
have been explored in this manner. 

Detailed investigations into mechanisms . - There have been several experiments 
directed to more specific and detailed information than available from the heated-tube 
experiments (refs. 66, 70, 73, and 113 and a paper entitled "Turbulent Heat Transfer 
to Carbon Dioxide Near the Critical Point" by R. J. Hanold, submitted to the 1970 
International Heat Transfer Conference). Probably the most useful details concerning 
the flow are velocity and temperature profiles. The only experiments in which profiles 
have been measured were performed by Wood (refs. 66 and 115) and by Wilson (ref. 113). 
In both cases, they surveyed radially across a vertically oriented heated tube near the 
exit. Wood used carbon dioxide at an L/d of 30. 7, and Wilson used hydrogen at an 
L/d of 128. 5. Wilson's data included hot wire measurements, in addition to tempera- 
ture and dynamic head profiles. In both cases, the major result was the appearance of 
the so-called M- shaped velocity profiles when the bulk temperature was near the trans- 
posed critical temperature, as shown typically in figure 21 and predicted analytically by 
Hsu and Smith. These are similar to the profiles shown in reference 116. Recent 
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Figure 21. -Typical M-shaped profiles as found by Wood and by 
Wilson (refs. 66 and 113). 


experiments of Bourke et al. (ref. 117) indicate the hot wire to be a useful tool in this 
regime. 

Kahn (ref. 73) and Hanold (in the paper mentioned on p. 40), both using carbon 
dioxide, conducted experiments designed to isolate some of the many interacting vari- 
ables and to highlight the influence of a single parameter. Kahn's method (refs. 73 
and 118) was to flow carbon dioxide at its critical point between parallel plates of large 
aspect ratio, with one plate heated and the other cooled at exactly the same rate so that 
no increase in the heat content of the fluid occurred. The heat-transfer coefficient 
peaked sharply when the hot wall (which was always the upper wall) approached the trans- 
posed critical temperature. The peaks became sharper as the temperature difference 
between the walls was decreased. Hanold used a horizontal rectangular channel with one 
wall heated. By changing the heated wall from upward to downward, he could compare 
transverse body-force effect. He found a negligible difference between the two cases 
and attributes it to a dominance of the flow by Kelvin- Helmholtz instability as opposed to 
Raleigh-Taylor type. 

A different result was observed by Hauptmann (refs. 70 and 119), who flowed near- 
critical carbon dioxide over a horizontal heated flat plate. Among his results are an in- 
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dication of a 30-percent increase in heat transfer for the same conditions with the heater 
oriented upward. Hauptmann's major contribution was a set of excellent, color, high- 
speed motion pictures coupled with good heat-transfer data. His data have already been 
discussed in figure 10. 

Shiralkar (ref. 74) and Jackson (ref. 120) have reported wall temperature spikes 
with near-critical carbon dioxide, which are discussed in the next section. 

In summary of the current, forced convection work, the following general observa- 
tions can be made: 

(1) Proximity to the critical point produces a sharp response in wall temperature, 
but the nature of the response is not uniform and the question of increased or decreased 
heat transfer is unclear. 

(2) There is evidence of a strong similarity to boiling and there seems to be little 
discontinuity in going from high-pressure boiling to the supercritical region. However, 
since the fluid is not two phase, the formulation is somewhat arbitrary. 

(3) Analyses and experiments directed specifically at critical- point heat transfer 
have produced some very interesting variable groupings which deserve continued in- 
vestigation. Among these are 

(a) Average specific heat, and integrated average properties, for example, 



(b) Minimum Prandtl number, (P r ) m i n 

(c) Eckert reference temperature parameter 



(d) Turbulence enhancement factors (see the section Theoretical Considerations 
in Force Convection) 


d(lnp) 


d(ln u ) 

dy + 
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* 

(e) Pseudoboiling parameters such as Stermann number and Martinelli param- 

eter 

(f) Stability limit (wall temperature spikes) defined by q = f(G) (see fig. 20) 

(4) Despite the fact that analysts developing the various correlation or mixing length 

methods have been able to substantiate their work with data, subsequent investigators 
have usually found it necessary to modify the method to fit their own data. No analysis 
is currently widely accepted nor has any analysis been thoroughly tested against the 
available data. 


Near-Critical Heat Transfer in Relation to Conventional Geometric Effects 


Curved tubes . - Systematic studies of the effects of curvature on a near-critical 
fluid are lacking. W. S. Miller (Rocketdyne Division, North American-Rockwell) has 
recently completed such a study, but the data are not yet available. Several studies have 
been made at higher pressures and for gases. The basic effect of curvature is to in- 
stigate a secondary flow such that the boundary layer is thinned at the concave surface 
(see fig. 22), and thickened at the convex surface. The flow patterns indicate that the 
core moves toward the concave surface and then flows back along the periphery to dump 
in the region of the convex surface. While the classical works of Ito (ref. 121) for tur- 
bulent flows and of Dean (ref. 122) for laminar flows deal with single-phase fully de- 
veloped fluids, the results of their analyses have been applied with varying degrees of 
success to the near-critical fluids (refs. 123 to 126 and a personal communication from 
W. S. Miller). 

Basically, Ito found that, for high- Reynolds -number fluids, the average friction 
factor for a curved tube was increased over that for a straight tube by the relation 


cr 

st 



(27) 


Hendricks and Simon (ref. 123) pointed out some effects of curvature on near- 
critical hydrogen, which are summarized as follows: 

(1) Rather conventionally, at high Reynolds numbers the concave surface enhances, 
and the convex surface degrades, heat transfer. The magnitude of these effects de- 
pends on fluid conditions, curvature-to-tube- radius ratio, and angular position along 
the bend. 

(2) Visual studies with small-bend- angle tubes, using liquid nitrogen, indicate this 
fluid is centrifuged to the concave wall rather than swirled in the normal secondary flow 
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Assumed velocity profile— v 




(b) Two-phase or near-critical flow. 


Figure 22. - Models of fluid flow in curved tubes. (Information on liquid film boundary 
received from K. J. Bell, Oklahoma State University.) 


patterns. Visual studies with steam indicated that droplets impinge on the concave wall, 
flow around to the convex surface, and are carried downstream. 

(3) Entrance conditions, profile similarity, and fluid history greatly influence the 
heat-transfer coefficients. Furthermore, the effects of curvature persist downstream 
of the bend and appear to be propagated upstream as well. 

The magnitude of these effects appears to depend on whether the wall is ’’wetted. ” 
For example, little circumferential variation in h was found for some subcritical pres- 
sure hydrogen tests at low T /T^. 

In subsequent tests at supercritical pressure and high T w /T b , it was found that h 
could vary as though the fluid was oscillating around the tube in some harmonic manner. 
Such oscillations were apparently instigated by a nonuniform heat-flux pattern at the 
upstream heating flange. They could also be induced by other factors, such as a non- 
uniform inlet velocity profile, or by the secondary flow itself as it moves from the core 
region along the periphery back to the core region. 

The effects of radius of curvature to tube radius R/r for hydrogen are shown in 
figure 23. At large R/r, the effect is small and approaches that of the straight tube. 
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log R/r 


Figure 23. - Concave-surface heat -transfer coefficient as function of ratio of radius of 
curvature to tube radius. 


At smaller R/r, the effects are quite large (2:1), but data seem to indicate two possible 
paths at lower R/r, as shown in figure 23. Which one is correct awaits experimental 
verification. 

McCarthy et al. (ref. 75) investigated near-critical nitrogen tetroxide and noted a 

nonuniform surface temperature phenomenon along the convex surface and a persistence 

of the curvature effects in the downstream region. Their experiment was complicated 

by the dissociation of nitrogen tetroxide; however, they still found up to a 2:1 increase 

in h due to curvature. In the work reported in reference 124, Thompson and Geery, 
cv 

investigating liquid hydrogen at P/P c ~ 5, also found a 2:1 increase in heat transfer 
because of curvature; however, they found no effects due to asymmetrical heating on 
noncircular flow passages. 

McCarthy, et al. (ref. 125), using hydrogen, found the heat-transfer coefficient to 
increase as 

h cr“ h st( 1+0 - 1 f) J £l ° < 28 > 


up to (L/d) = 10 and to become fully developed thereafter. He suggested the following 
equations: 
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where (Nu) g j. is given by table IV-C of reference 55, and by equation (25). 

Taylor (ref. 126) investigated the data of references 124 and 125 and found they 
could be adequately correlated by using Ito's parameter: 


( Nu W _ jO. 05 
< N “>calc ’ b 


Concave 


(31) 


( Nu )cx j-0. 05 
< Nu >calc ’ b 


Convex 


(32) 


where (Nu) ca ^ c is given as equation (9). While equations (28) to (32) were developed for 
hydrogen, similar correlations could be anticipated for other fluids, at high T /T^ 
values ("nonwetted" walls). "Wetted” walls and boiling have been investigated by 
K. J. Bell of Okalahoma State University. He found a substantial increase in h on the 
convex wall for wetting two-phase fluids at low T (private communication). The 
visual steam flow patterns noted by Hendricks and Simon lend credence to these results. 

Twisted tapes and rifle boring . - If the bulk fluid is near the wall, an augmentation 
in the heat-transfer coefficient will occur. Bartlit and Williamson (ref. 127) induced 
swirl flow in liquid hydrogen using a twisted tape in a long tube and found the heat- 
transfer coefficients to be approximately those predicted by the standard Dittus-Boelter 
equation. This is in contrast to the same pipe without a swirl inducer, where the ex- 
perimental h was much less than the predicted h (fig. 24, see also ref. 128). 

E. Michaud and C. P. Welsh (Badcock and Wilcox, Barberton, Ohio) have demon- 
strated a similar effect with near-critical water (personal communication, see also 
ref. 129). The lands of a rifle-bored tube both induce swirl and trip the flow; these 
geometry effects eliminated the wall temperature "spike" noted in a similar smooth 
tube under the same conditions (fig. 25). 

The enhancements, of course, require a pressure-drop penalty. J. H. Van Sant 
(Laurence Radiation Laboratory, University of California, Livermore, Calif. ) analyzed 
the results of Kidd (ref. 130) and found that the enhancement in heat transfer was at the 
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direct expense of pressure drop (personal communication). The optimization of heat 
transfer, pressure drop, and heat- exchanger compactness requires much effort. 

Body-force orientation . - Orientation of the heat-transfer apparatus with respect to 
gravity can have considerable influence on the results. Hall, et al. (refs. 131 and 132), 
Jackson and Evans- Lutterodt (ref. 120), Wood (ref. 115), and Shiralkar and Griffith 
(ref. 74 and 128) have found that parallel buoyancy effects can substantially alter the 
velocity profiles. Shiralkar et al. (ref. 74) and Jackson and Evans- Lutterodt (ref. 120) 
found wall temperature spikes in carbon dioxide at elevated fluxes for upward flow. In 
reference 120 the spikes were eliminated by downward flow. The effects also seem 
quite dependent on whether the bulk temperature is above or below T*. However, the 
present authors found no noticeable effect with high-velocity, near- critical hydrogen in 
up- and downflow configurations. The first-order answer lies in the ratio of buoyancy 
to inertia forces which, for the authors case, was quite small (Gr) f /(ReL = 4x10 . 

1 I O _ Q 

Shiralkar and Griffith (ref. 74) found the effect to be noticeable when Gr/(Re) > 4x10 
" 2 

for water and 10 for carbon dioxide. These values and conclusions are in agreement 
with the computations of Hsu and Smith (ref. 133). 

Transverse effects, studied by Hauptmann (ref. 70) and Hanold (in the paper cited 
on p. 40), have already been discussed. 

Entrance effects . - Papell and Brown (ref. 134) investigated entrance effects in 
near-critical hydrogen. The experimental h at high heat fluxes were 1. 5 times higher 
for an undeveloped entrance than for a hydrodynamically developed entrance, as illus- 
trated by figure 26. This effect persisted an abnormally long distance down the tube. 
Little effect was found for the low-heat-flux case. 



Figure 26. -Comparison of developing flow to fully developed flow for heat transfer in the en- 
trance region with near-critical hydrogen. (From ref. 134.) 
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In summary, the geometric effects of significant interest in near-critical heat trans- 
fer are 

(1) An increase of about 2:1 in the heat-transfer coefficient on the concave side can 
be expected because of curvature for gases and fluid systems at high T w /T b * This en- 
hancement seems to be a centrifuging effect, particularly in small-bend-angle tubes. 

(2) The magnitude of enhancement (or degradation on the convex side) appears to be 
related to the Ito parameter. As a reference for computing enhancement, for the "dry- 
wall’' case, correlations of the type given in equation (9) appear to be the most satis- 
factory. For the "wet-wall" case, the augmentation can shift to the convex wall. 

(3) Devices which centrifuge the high-density core to the wall or break up the low- 
density boundary layer, such as twisted tapes and spiral grooves, enhance near-critical 
heat transfer significantly. 

(4) Experiment and analysis suggest that a parallel body force should begin to in- 

2 -3 

fluence heat transfer when Gr/(Re) >5x10 

(5) Tube entrance effects persist much further down the tube in near-critical fluids 
than in gases. 


Theoretical Considerations in Forced Convection 

In the free and natural convection areas, reasonable success has been realized by 
using conventional variable- property approaches. For example, in natural circulation 
loops, the arguments of Schmidt, et al. (ref. 6) have stood the test of time and can be 
seen in use today in such correlating equations as equations (17) and (18). In forced 
convection, the fluid is not allowed to naturally seek its own circulation patterns. Con- 
sequently, the complex flow interactions caused by the rapid expansion of the near- 
critical fluid are superimposed on the normal forced convection patterns, and conven- 
tional forced convection approaches seem inadequate. A fair amount of analytic atten- 
tion has been turned toward explaining these phenomena. This will be briefly recapped 
in this section. 

Mixing length analyses . - By far the most widely used approach is the Prandtl 
mixing length concept. The basic idea is that the Reynolds stress terms of the Navier- 
Stokes equations can be made to look like laminar shear terms by the introduction of a 
turbulent viscosity or eddy diffusivity. The shear and heat flux expressions assume the 
popular form: 


T 


du 

= V — +pe 


dy 


du 

dy 


(33) 
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(34) 


q = 



' h 


dT 

dy 


Subject to the following nondimensionalization, Deissler (refs, 
these equations apply throughout a supercritical flow field: 



y + = n!y 

V 

o 


T + 




0 = 




where u* is the well-known shear velocity 



The nondimensional equations become, for = e 


2- = flL + A 

T o \^o p o ) dy + 


q_ 


X 

X o( Pr )o 


P 0 



dT + 

dy + 


135 and 136) proposed that 


(35) 


(36) 

(37) 
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Deissler (ref. 135) numerically integrated these equations for water at P/P = 1. 56, 

c 

with moderate success. He suggested using a reference temperature T x (discussed in 
the section Heat Transfer in Forced Convection Systems) for evaluating properties, and 
Eckert (ref. 99) consolidated the data spread by introducing the following dimensionless 
temperature: 


E = 


T *- T b 



(38) 


In that same year, Goldman (ref. 137) introduced a technique to transform the 
variable- density problem to the constant- density form. This was done by using a vari- 
able density in the shear velocity: 


(39) 

The form of equations (36) and (37) remains unchanged, but all the nondimensional vari- 
ables are changed. They are usually designated u ++ , y ++ , and T ++ to indicate the 
change (see table V). Goldman also suggested a different form for the heat-transfer 
coefficient, which evolves from the Dittus-Boelter equation (ref. 138): 

n H 0. 2 

h' = — (40) 

G°.S 

These modifications have been popular, but universal success has not been achieved. 

Most of the research in mixing length theories for near-critical fluids has centered 
on adapting the eddy diffusivity to comprehend the near-critical phenomenon. To review 
the basic eddy diffusivity concept, look at two prominent examples, those of Van Driest 
(ref. 139) and Deissler (ref. 135): 

Most eddy diffusivities e can be placed into two categories: (1) continuous, and 
(2) multiple part. In the continuous case, Van Driest (ref. 139) argues that one form of 
e should apply over the entire regime; while for the multiple-part case, Deissler 
(ref. 135) points to a near-the-wall region and an away-from-the-wall region. These 
and other eddy diffusivities are present in table V. 

Van Driest continuous approach: 



(41) 


£ = K 2 y 2 [l - e-^f 


du 

dy 
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TABLE V. - EDDY DIFFUSIVITY EXPRESSIONS 


Dimensionless eddy diffusivity, 


;_ c -n 2 ( V ">uVj nV 


Constraint Energy-momentum Reference 

solution 


e, r - 2 Coupled Deissler and Presler (ref. 213) 

im compressible n -0.124 


-K 2 ( V T) (du + dy + ) 3 /(d 2 u + dy +2 ) 2 l ? 


1 - e -<y + ^K 2 48*^ 


u? re u + - f pL du + 

j o y "o 

y ^ y -.-/ y+ 

J 0 \"o ” 


du _ _*o du 
dy + r] dy -1 " 


JL - pL (l2 

,* 2 


Decoupled Van Driest (ref. 139) 

incompressible K = 0. 4 A* = 26 

Coupled transform Goldman (ref. 78) 
of compressible 
to incompres- 
sible 


(Note: t y D 


-0.0285 v ■ v. 


pj where A = 30.2e 

d In p 

(1 + F ) = e + F - d y + 


Hess and Kunz (ref. 110) 


Hsu and Smith (ref. 133) 


1 + F h> ■ ‘J F h 


; 1 m, 

-11 tanh V-} 
llJ 


r jz« + 

l/o lfPo J 

7y + fl - o' 0 002 *y t ) 2 | 


Decoupled 


Hall, Jackson, and Kahn 
(ref. 118); r(y) is linear; 

A = 0.4; du ++ /dy ++ ; conven- 
tional velocity profile 

Petukhov and Popov (ref. 214) 


Popov (ref. 215) 
Reichart (ref. 216) 


Melik- Pashaev (ref. 217) 


Tanaka, Nishiwaki, and 
Hirata (ref. 65) 



Deissler multiple-part approach: 


£ = n 2 uy^l - e-" 2 ^ £ < y + * 2 




e = K 


du 

dy 


f djj' 

Wy 5 


1 _ e -(K 2 /^)(du/dy) 3 /(d 2 u/dy 2 ) 


( 42 ) 


y + ~ 2 


Equations (41) and (42) were developed primarily for constant-density analyses. They 
have been used successfully for variable-property gases (Deissler, ref. 135). The near- 
critical research has required modification of these basic forms. These modifications 
are summarized in table V. A few are discussed in this section. 

Hess and Kunz (ref. 110) found that, for near-critical hydrogen, the damping factor 
of equation (41) had to include a term dependent on kinematic viscosity. Their work in 
this area led to the correlating equation most frequently used for near-critical hydrogen 
(eq. (25)). 

Hsu and Smith (ref. 133) proposed that, in the presence of a strong density gradient, 
shear should be represented as the rate of change of momentum. The resulting eddy 
diffusivity expressions may be written in terms of equations (41) and (42) as 


£ m = < l + V £ 




where e is given by equation (42), and 

d In p 

f - d y + 

m d In u + 

dy + 


F h = 


d In p 
dy + 

d In CY 1 


dy 4 


. J 


(43) 


(44) 


Equation (43) reduces properly for constant density, and either form of e (eq. (41) 
or (42)) may be used. Hsu and Smith (ref. 133) also performed a first-order force 
balance on a differential volume to derive an equation relating shear, buoyancy, and 
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pressure drop. The success of their equation in correlating the data of Bringer and 
Smith (ref. 100) indicates that such a modification is warranted, and is another step in 
the proper direction. The main problem here seems to be that, while each researcher 
has success correlating his own data, subsequent researchers had troubles. Further, 
Hall (ref. 131) using several different techniques on the data of Koppel (ref. 140) also 
found a degree of insensitivity in the theories, as illustrated in figure 27. 

Hall (ref. 132) developed a model for flow in a vertically oriented pipe or channel 
in order to demonstrate some of the buoyancy effects which have been found with near- 
critical fluids (see the section Body-force orientation ). Halls' theoretical approach, 
while limited and analogous to the two-phase or two-component flow approach, merits 
attention because it isolates some governing parametric forms. The first model as- 
sumes a discontinuous change in transport properties and density at y = 6, with ''gas- 
like” properties at y < 6 and liquid properties for y > 6 (e. g. , see fig. 28). The 
shear distribution is a two-part linear form (see fig. 28). The velocity profile is lami- 


Exptrimtnt 

Correlation 

Thtory 


o D«iss!«r 11954) 
a Goldmann (1954) 

□ W«idtracht and Sonn«mann (1960) 

• Miropolsky and Shdsman (1957) 

■ Petukhov. Krajnoschekhov and Protopopov (1961J 

* Kutateladzt and Itontiev (1964) 



Figure 27. - Comparison of various correlations and 
mixing length theories to carbon dioxide data of 
Koppel. Ratio of pressure to critical pressure, 
1.028; mass flow rate, 92. 2 grams per centimeter 
per second; heat flux, 22. 5 joules per centimeter 
squared per second; characteristic diameter, 
0.493 centimeter. (Reproduced from ref. 131.) 



(b) Shear distribution. 


Figure 28. - Property variation and shear distribution for the models of 
Hall (ref. 132). 


54 



nar for y < 6 and follows the conventional velocity defect law for y > 6. With these 
assumptions, the momentum equation 


dz _ Ajo 
dy AX 


PS 


where Ap/AX = Constant, can be reduced to 



w 





Here b„ = y + evaluated at the tube radius or channel half width, £ = 5/b„, and 
11 3/2 H 

(Gr)' = P w (P core “ Pw^hAw is conventional Grashof number as defined for two- 

phase or two-component flows. 

In the second model, Hall (ref. 132), employed equations (36) and (37) with 
Van Driest's eddy diffusivity, equation (41). However, the shear stress distribution 
t/t q was assumed to be of two parts (see the Deissler approach, eq. (42)): 


-H=1-JL 


T 


O 


H 



near the wall 



in the core 


where A, = t./t . Hall (ref. 132) found little difference in the results given from the 

A o w 

two models, which reinforces an earlier remark regarding the insensitivity in the 
theories. Hall and Jackson (ref. 141) and Leontiev (ref. 142) have tried to relate this 
phenomenon to that of lamarization. 

Each theoretical technique has a great deal of merit and correlates the data pre- 
sented in each author's respective papers. Using near-critical hydrogen, we at NASA 
and also Szetela (ref. 143) have tried the techniques of Deissler (refs. 135 and 136) and 
Goldman (ref. 137) with very limited success. The technique of Hess and Kunz (ref. 110) 
was tried, and success was limited, in the main, to those data with bulk temperatures 
above the T* regime. The integral average property approach of Goldman (ref. 137) 
gave similar limited results. The modifications of Hsu and Smith (ref. 133) and Brokaw 
(ref. 109) were also tried again, with only limited success. Modifications of the modi- 
fications, such as using different dampening factors, including upstream effects, etc. , 
were also attempted; however, the results were very disappointing. 
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The solution of equations (36) and (37) (or modifications, eqs. (43) and (44)) involves 
the iteration of a numerical integration to attain a solution. The iteration is required 
to satisfy the continuity equation and to determine one parameter. In pipe flow, this 
parameter can be 0 = r 0 T 0 > V or r+ * However > r+ is usually fixed at the 

pipe centerline and t q is considered to be known; thus q Q and T q remain as adjust- 
able constants, within the parameter 0. In the boundary- layer flows, at the edge of the 
boundary layer 5 + , either 0 or r Q can be considered to be the parameters; usually 
<5 + is altered at a fixed j 3 and t q to give the solution. The point is that for near- 
adiabatic flow of a gas there exists an adjustable parameter to force a solution to the 
problem. However, this freedom becomes compounded when applying these equations in 
the near-critical region because the following factors are not known: 

(1) Does e/e h — 1? (Some of the best results for near-critical fluids in a pipe 
attained by the authors (unpublished computer programs by Robert C. Hendricks) were 
for e/e h > 1. 


Curve Theory of- 

1 Oeissier (ref. 135) 

2 Goldman (ref. 137) 


Empericai equations of- 

3 Miropulskii and Shitsman 

(ref. 103) 

— - — 4 Krasnoshchekov and Protopopov 

(ref. 145) 

— ^ — 5 Swenson, et. al (ref. 61) 


30xl0 3 



Figure 29. - Heat-transfer coefficient as 
function of wall temperature for water. 
Pressure, 345 atmospheres (at 417° C); 
ratio of wall temperature to film tem- 
perature, 1.25; mass flow, 2150 kilo- 
grams per square meter per second; 
characteristic length, 9.4 millimeters. 
(From ref. 144.) 
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(2) Are the expressions for e or correct? 

(3) Is the functional form of t/t q and q/q Q realistic (i. e. , constant, linear, 
cubic, etc.)? 

(4) Do the convective terms in the momentum and energy equations significantly 
alter the results? These terms are considered in the next section. 


In view of these unknowns, several different velocity and energy distributions may in- 
tegrate out to give the same solution. However, the convergence may be quite slow and 
difficult to achieve; and furthermore, the result may be meaningless. 

Some of these difficulties have also been encountered by Petukhov. In a recent sur- 
vey article (ref. 144), he lists three basic difficulties in attaining solutions for heat 
transfer in the near-critical regime: 

(1) The nonlinearity of the governing equations with variable properties 

(2) Ignorance of the laws of turbulent heat transfer with variable physical prop- 

erties: To determine these laws will require detailed profile measurements in 
the near-critical regime - a task of considerable difficulty. 

(3) Difficulty of interpretation and generalization of near-critical data because heat 

transfer and pressure losses are functions of many variables 
In figure 29, Petukhov (ref. 144) compares the theoretical results of Deissler (ref. 135) 
and Goldman (ref. 137) and the empirical equations of Miropulskii and Shitsman 
(ref. 103), Krasnoshchekov and Protopopov (ref. 145), and Swenson, et al. (ref. 61). 

Petukhov (ref. 44) points out that the theoretical analyses are in considerable dis- 
agreement with the experimental results although no data are given in figure 29 for com- 
parison. In reference 144, Petukhov also presents an analysis of Popov which is valid, 
unfortunately, only for carbon dioxide data. Petukhov (ref. 144) classifies the near- 
critical heat-transfer regimes as 

(1) Normal: where the dependence of Nu on Re and Pr is approximately the same 

as for the constant properties case 

(2) Reduced heat transfer: where T buUf < T*, and certain combinations of mass 

flow and heat flux lead to a sharp decrease in heat transfer and a rapid in- 
crease in the wall temperature 

(3) Improved heat transfer: where T bulk < T*, with oscillations of sufficient 

amplitude ("whistling”) and high frequency (1400 to 2200 Hz) at higher heat 
fluxes to cause a flattening of the axial wall temperature profile 
These classifications are in agreement with those presented herein but, perhaps, are 
formulated in a different manner. 

Acceleration - strain rates. - Deissler (ref. 146) found the strain rate (acceleration 
or deceleration) to be important in turbulent heat transfer. For example, the strain ef- 
fects of longitudinal acceleration on homogeneous turbulence intensities are shown in 
figure 30. The ordinate is the square of the turbulence intensity of the transverse veloc- 
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Axial position, X + ■ X/( -Xg) 

Figure 30. - Effect of acceleration or 
longitudinal strain on heat transfer 
for positive strain rates. (From 
ref. 146.) 


ity component normalized with respect to an initial intensity condition (v')V u ^J/ 

[ (v') 2 /i?|. It can be shown that this dimensionless parameter is related to the square 
u 0 ' Uj 


of a Stanton number ratio (St)/ (St) Q where the denominator is the initial Stanton number 
before acceleration 
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(v ) 2 


0 


u. 



(45) 


The abscissa of figure 30 is axial position; the parameter s + is a strain rate. Fig- 
ure 30 shows that the local intensity, or Stanton number, progressively decreases when 
subject to a strain rate. 

Under heating conditions, a near-critical fluid would undergo large longitudinal ac- 
celerations while expanding in a constant-area pipe. Experiments by Weiland (ref. 147) 
and Taylor (ref. 148) with gases at high heating rates suggest that such expansion effects 
occur. Deissler (ref. 146) showed that the data trends in references 147 and 148 could 
be explained by strain or acceleration effects. Such an analogy could also be drawn for 
the laminarization data of Moretti and Kays (ref. 149). Perhaps some of the unusual 
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Figure 31. - Concept of simultaneous hydrodynamic and penetration mechanisms. 
(From ref. 150). 


distributions of heat-transfer coefficient noted for near-critical hydrogen (refs. 41 and 
134) relate to strain effects. 

Penetration model. - An analogy between heat transfer in fluidized beds and near- 
critical forced convection has been suggested by Graham (ref. 150). Using the penetration 
theories of fluidized beds, Graham proposed a penetration model for near-critical fluids. 
The mechanism (see fig. 31) was assumed to be a disruption of the sublayer by the peri- 
odic migration of fluid packets from the outer edge of the boundary layer to the wall, 
where they absorb energy in a transient manner. The average heat-transfer coefficient 
can then be expressed as a linear sum of the conventional film-heat-transfer coefficient 
and a penetration coefficient with the proper area weight fraction assigned to each 



(45a) 


where 
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(45b) 



contact 


The theory suffers from the lack of a theoretical means of determining T con t ac t> 
but has the distinct advantage of reducing to the proper expression outside the near- 
critical region and integrates the physics of the near- critical problem in a realistic 
manner. The technique has been used to correlate near-critical hydrogen, carbon di- 
oxide, and water data. 

In summary, the theoretical investigation in forced convection in the near-critical 
region suggests 

(1) Conventional mixing length theories must be modified to account for the tremen- 
dous fluid expansion in the near-critical region. These modifications could be in the 
equations, following Goldman's example (ref. 137); or in the eddy diffusivity, following 
Hsu and Smith (ref. 133); or in forming a penetration mechanism, Gral am (ref. 150). 

(2) Applications of these and other modifications of mixing length theory to near- 
critical data have met with limited success. The theory will probably also have to be 
made to comprehend upstream history, acceleration, etc. , in order to achieve wide 
success. 

(3) Further analysis of accelerated flows should consider the effects of strain rates. 
While the theory of Deissler predicts the proper trends, the magnitudes associated with 
a real near-critical fluid are not correct. Thus, the boundaries and the variation in 
fluid properties must be properly considered (not an easy task). 


Oscillations 

General remarks . - Fluid oscillations play a major part in near-critical phenom- 
ena, and several heat-transfer studies (refs. 42, 75, 98, 125, and 151 to 153) have re- 
ported some type of instability in natural and forced convection systems. The frequency 
of reported oscillations varies from 1/20 to 10 000 hertz over a considerable range in 
amplitude (as high as 0. 3 of test pressure). The higher- frequency oscillations, 300 to 
10 000 hertz, seem to be associated with longitudinal and lateral oscillations - a gas- 
like fluid behavior; and the lower-frequency modes are most common to the near-critical 
region. Thus, discussion is restricted mainly to the low-frequency spectrum. 

Fluid behavior in the two-phase and near-critical regimes appears to be similar, 
and different modes of oscillations have been identified and analyzed. The analyses have 
been successful in predicting the onset and frequency of oscillation but less so in pre- 
dicting the amplitudes. Flow excursions, as well as periodic flows, have been studied 
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Ratio of heat flux per unit area 
to single-tube liquid-oxygen flow 
rate, qA/G, kj/im^ikg) 


Figure 32. - Percentage of heat-exchanger-test data samples showing 
steady operation against heat flux per unit area for unit mass flow rate. 
Single-tube data coil pressure, 40 to 47. 5 atmospheres. (From ref. 39. ) 


by Zuber (ref. 39). Similar oscillations in natural and forced convection have been re- 
ported. Several authors after observing oscillations (refs. 39 and 41) have delineated 
correctable and noncorrelatable data based on mass flow rate and heat flux, as illus- 
trated in figure 32 (see also fig. 20). 

Oscillations in the near-critical regime have the common characteristic of persist- 
ence for a significant period of time. This means they occur at, or very near, a natural 
frequency of the system, and the damping forces are insufficient to quickly arrest them. 

Oscillations may be divided into four types, of which three may be separated by their 
frequencies. The other type is associated with the oscillation of a near-critical fluid be- 
tween laminar and turbulent flow. Due to the sensitive nature of near-critical fluids to 
fluctuations, the laminar-turbulent oscillations may be characteristic of any heat- 
transfer system of equal sensitivity. These oscillations were discussed in the section 
Heat Transfer in Free Convection Systems. 

The other three types of oscillations may be classified as (1) thermal-acoustic 
(5 to 300 Hz); (2) system oscillations (1/20 to 5 Hz), which are concerned with the be- 
havior of all the fluid in a loop or system; and (3) external oscillations (e. g. , in pipes, 
in fluid conduits, and in valves). This latter type would be induced and sustained by the 
equipment and have a frequency which is that of the static system. They will not be dis- 
cussed further. 

Examples of both thermal-acoustic and system oscillations are shown in figure 33. 
These oscillations appear to have a significant density change as the source of their be- 
havior, which is characteristic of the onset of the two-phase and near-critical regimes. 
Here, the density change becomes very sensitive to further additions in energy or 
changes of temperature. Generally, oscillations will occur when the bulk of the fluid is 


61 




TIME, SECONDS 

(a) Thermal -acoustic oscillations (ref. 40). 
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(b) System oscillations (ref. 151). 


Figure 33. - Examples of thermal -acoustic and system oscillations reproduced from references 40 and 151. 


in a liquefied condition and the boundary layer appears to be occupied by a lower-density 
fluid. Most analyses of oscillations assume density to change abruptly at the saturated 
liquid and vapor points and at the equivalent points in the near-critical region. This as- 
sumption tends to more sharply define the oscillating behavior. 

Thermal-acoustic oscillations . - The thermal-acoustic oscillations (5 to 300 Hz) 
may be assumed to occur from localized events, in contrast to those of the fluid-body 
type, which involve the motion of the entire fluid system. Thurston (ref. 152) identified 
two types of these oscillations in his forced convection system. The first type is a 
Helmholtz oscillation, which is associated with a cavity connected to a large plenum, 
such as the atmosphere, with a substantial restriction between the two elements. The 
natural frequency of this system may be expressed by the equation 


W n 



(46) 


This Helmoholtz resonance has not been reported in most of the systems concerned with 
this type of study and may be presumed to occur only when the system obviously fits the 
analytical assumption rather closely. 

The second type of acoustical oscillation is a sonic-wave resonance, identified by 
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Thurston and others, and is associated with the velocity of sound and the length of fluid 
in the flow loop or fluid system considered. The natural frequency of this system may 
be expressed by the following equation: 


W n = 


c 

2L' 


(47) 


The period for these oscillations is that time required for the acoustic wave to tran- 
verse one-half the fluid length. The frequency of these oscillations has been found ex- 
perimentally to be a very weak function of heating rate, and appropriately it changed 
with fluid length. 

Several authors have proposed mechanisms which would produce and reinforce these 
acoustic oscillations. These mechanisms examine the behavior of the boundary layer 
during the passing of a pressure wave. Thurston (ref. 152) suggested that such a wave 
would tend to compress the boundary layer, improve its thermal conductivity, and cause 
an increased heat-transfer rate. Immediately following this process there would be a 
tendency for the boundary layer to expand, with a subsequent decrease in thermal con- 
ductivity and decreased heat-transfer rate. This mechanism would reinforce the acous- 
tic wave as it proceeds through the field. 

Thurston, et al. (ref. 40) developed two similarity parameters to determine the 
onset and amplitude of the thermal- acoustic oscillations. One is a function of the 
Stermann parameter (boiling number) and the other a function of properties, fluid con- 
ditions, flow rate, and geometry of the system. 

Thurston, et al. (ref. 40) proposed that the onset of these oscillations can be em- 
pirically determined as a relation between the Stermann parameter and the specific 
volume number described by the following expression: 


Sr < 0. 0045(Svf 0,75 


where 


Sr 



and 



(48) 


(48a) 


(48b) 
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The locus of the onset of oscillations is described by the equality in equation (48) (Sr) Q 
and the stable points by the inequality. Other authors have identified this threshold by 
applying stability criteria to the fluid system, employing the general conservation equa- 
tions. 

The amplitude parameter developed by Thurston was roughly correlated as a func- 
tion of Stermann number 

Amplitude parameter: 


Am 


p*W n d 

Gh* 

fg 


(48c) 


Amplitude correlation: 


Am = 0. Olljsr - (Sr)J 0 ' 66 (49) 

System oscillations . - Thermal-ascoustic oscillations are confined to the heater 
section, while system oscillations involve the motion of the fluid in the entire loop. 
System oscillations (1/20 to 5 Hz) seem to be present to some degree in nearly all near- 
critical heat-transfer experiments. 

Thus, while the thermal-acoustic oscillations may have a significant amplitude and 
this, in turn, may affect the heat-transfer and general system behavior, most facilities 
will not fail to function because of their presence. System oscillations however, could 
very conceivably prevent the operation of the loop in several different ways (e. g. , a 
pressure variation effecting cyclic heat-transfer conditions to produce forces which may 
cause mechanical failure). The frequency of these oscillations is somewhat lower than 
for the thermal-acoustic oscillations. Also, they are one function of the energy input 
to the fluid, and another function of the system geometry. 

Analysis of system oscillations, whether in a two-phase or in a supercritical re- 
gion, tend to assume a boiling-like phenomenon with substantial density changes. The 
basic approach is to supplement the time- dependent, one-dimensional conservation 
equations (continuity, energy, and momentum) by an equation of state which exhibits a 
discontinuity at the saturated liquid line or at the equivalent point in the near-critical 
regime. Two general methods have been used to predict the onset of these oscillations. 
The first, attributed to Craya and Boure (ref. 154), is called the density effect model. 
This method has been used and extended by a number of authors (refs. 155 to 157). 
Walker and Harden (ref. 156) found six independent variables which would influence the 
inception point. Three of these are concerned with geometry and the friction factor 
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which could be considered essentially constant for any given system. The remaining 
three variables are fluid velocity, gravitational force, and subcooling. The critical 
frequency, the same as that of an oscillating manometer, is given by 


1 

2g(p f - Pg ) 

2n 

c(p f + p g ) 


( 50 ) 


Walker and Harden (ref. 157) were able to compute a threshold surface involving the 
three variables and to predict instabilities in a Freon- 114 natural convection loop with 
good reliability. 

The second criteria used the Nyquist-type of test to determine the stability of the 
system. This method is widely used for electrical systems and can be employed in fluid 
systems by utilizing analogous behavior. For example, the resistance and capacitance 
of fluid system components can be modeled electrically. Examples of such modeling can 
be found in the work of Zuber (ref. 39), Friedly et al. (ref. 158), and Wallis and Heasley 
(ref. 159). Examples of working heat- exchanger oscillations are found in reference 160. 

The following remarks can be made in summary: 

(1) Oscillations in near-critical flow can be divided into four groups: laminar to 
turbulent, thermal acoustic, system, and external. 

(2) The large density changes of a near-critical fluid are the basic source of the 
oscillating behavior. 

(3) The critical frequency for the thermal- acoustic and system oscillations can be 
predicted by utilizing simple concepts of mechanics. 

(4) Conventional stability criteria may be employed to predict the thermal-acoustic 
and system oscillations with a fair degree of reliability. 

(5) Attempts to predict the amplitude of the oscillations has met with only limited 
success. 


Choking Phenomenon 

Although the authors are aware of only one reported work and their own current 
work on this subject, choking flow in the near-critical region may be unique and of con- 
siderable interest. Choking flow, often referred to as critical flow, may be defined as 
the condition where the mean velocity and the mass flow rate are at a maximum with re- 
spect to further reductions in the downstream pressure (refs. 161 and 162). 

W. G. Zinman, in a private communication, points to the observed similarities in 
various critical phenomena (e. g. , choked flow, detonations, shock waves, critical points 
of ferromagnets, thermodynamic critical point) and suggests there should be some com- 
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mon theoretical base for each system. For example, the thermodynamic critical region 
and choking flows are essentially unstable phenomena, and experimental measurements 
to determine either are destined to some degree of failure. The experimental problems 
encountered in obtaining data for either of these critical conditions are monumental. As 
critical conditions are approached, changes occur very rapidly. In the thermodynamic 
critical region, these are property changes; and in choking flow, they are usually char- 
acterized by pressure changes in the flow direction. Because of these very rapid 
changes, relaxation times become important. It follows that for many systems, at the 
rates of change required to maintain equilibrium, it is unlikely that equilibrium will be 
achieved nor can the system be assumed to be frozen. Therefore, the problem of 
choking flow near the thermodynamic critical point is compounded by the necessity of 
considering two critical and elusive phenomena combined. 

For a single-phase, homogeneous fluid, the critical flow rate may be described by 

) (51) 

process 

If the process is assumed to be isentropic, the choking velocity is equivalent to the sonic 
velocity. 

If it is assumed that the fluid is homogeneous in properties and velocity and in ther- 
mal equilibrium at the critical point, the choking and sonic velocities as described by 
u = dP/dp would be zero. This is because the derivative in equation (51) goes to zero 
at the critical point. Also, in the near-critical region, the value of this derivative tends 
to be quite low. Some experimental work reported by Williamson and Chase (ref. 163) 
gives the measured sonic velocity of helium in the near-critical region. This is shown 
in figure 34. Examination of the figure shows that, indeed, the velocity does appear to 
approach zero near the critical point and to be at a relatively low minimum in the region 
of the transposed critical. Van Dael, et al. (ref. 164) reported similar results for 
argon. Further, they indicate a practical lower limit on sonic velocity of about 
100 meters per second. 

It would appear that if, indeed, the very low choking velocities were probable in the 
near-critical region, some experimental measurements (a discontinuous pressure pro- 
file accompanied by abnormally large pressure loss, or visual observations) would have 
been reported which would indicate the presence of choking flow in the system. Since 
none of these conditions has been reported in the near-critical fluid flow and heat- 
transfer studies, it must be assumed that the choking velocities are somewhat higher 
than those indicated by the sonic velocity or by studies which assume thermal equilib- 
rium with a homogeneous fluid. The work of Hesson and Peck (ref. 165) with carbon 
dioxide appears to be the closest to the critical region. Hesson and Peck mapped two- 
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Figure 34. - Overall view of sound velocity of helium-4 as function of pressure near 
critical point. (From ref. 163.) 


phase critical flow for stagnation pressures up to the critical; however, no data were 
taken above critical pressure or for throat pressures in the near-critical region. A 
similar experimental program for nitrogen flowing in a converging-diverging nozzle is 
currently under way at Lewis Research Center (unpublished data by R. E. Henry, 

R. J. Simoneau, R. C. Hendricks, and R. Watterson). Preliminary data have been 
taken for stagnation temperatures ranging from 85 to 128 K and for stagnation pressures 
in the range of 26. 4 to 41. 8 atmospheres. These results indicate that the choking phe- 
nomenon in this region is characterized by a nonequilibium expansion. In addition to 
the conditions listed, the stagnation pressure range will be extended to 68 atmospheres 
in an attempt to create near-critical conditions at the nozzle throat. 

For the two-phase-flow case, which should represent an extreme of the near-critical 
two-species choking flow, it has been shown that assumptions of nonequilibrium or of 
different phase velocities are necessary to explain the experimental data (refs. 166 
to 168). 

In summary it can be said that 

(1) If a near-critical fluid were in thermodynamic equilibrium, the sonic velocity 
would be low and the flow might be expected to choke very easily. 

(2) None of the normal indications of choking have been observed in the near-critical 
experiments. This suggests the possibility of thermodynamic nonequilibrium in these 
flows, which is also consistent with two-phase choking experience. Preliminary nozzle 
data for near-critical nitrogen indicate the expansion to be nonequilibrium. 
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Zero-Gravity Operation 


The authors are aware of only one qualitative heat-transfer study of a near- critical 
fluid at very low gravity. Potter and Brill (ref. 169) placed 2. 5 liters of liquid nitrogen 
on an Atlas ballistic flight. They present data for nearly 1500 seconds, during which 
time the gravity force increased from 6x10"^ to 0. 1034 g’s. Liquid-nitrogen pressure 
level was controlled by an internal heater and periodic flows. Although Potter and Brill 
report that zero-gravity operation of the vessel was not significantly different from 
earthbound tests, their pressure profiles do indicate some differences. Data in this 
area are lacking and work should be performed. 


SUMMARY OF RESULTS 

Various areas of research on near-critical heat-transfer have been surveyed. The 
results for each area have been summarized at the end of that particular section. These 
summaries are repeated here for the convenience of the reader. 


Fluid Properties 

1. Classical Van der Waals models do not accurately describe the state of a near- 
critical fluid, but may be satisfactory in first-order considerations. 

2. The discrepancies in the law of corresponding states are accentuated in the crit- 
ical region, and the law should be used only between fluids of similar molecular struc- 
ture. 

3. A thermally disturbed near-critical fluid requires long relaxation times to return 
to equilibrium. 

4. Equations of state of the conventional virial type can be made to accurately map 
a near-critical pressure-density-temperature surface; however, they are not good for 
thermal properties requiring derivatives of the equation of state. 

5. A new analytic equation of state involving chemical potential, density, and tem- 
perature, based on the Ising lattice gas model, has been developed for the near-critical 
region. 

6. The evidence for the existence of a strong anomalous spike in thermal conduc- 
tivity at the critical point is overwhelming. The spike can be estimated by treating the 
fluid as a dissociating polymer and superimposing the results on the conventional 

X - X* = f .(p) curve (eq. (5)), where X is thermal conductivity, X* is the thermal con- 
ductivity at low pressure, and f ^(p) is a function describing variations in thermal con- 
ductivity. 
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7. If a similar spike occurs in viscosity, it is sufficiently weak to be ignored. 
Thus, 7] - T)* = fgfp) (eq. (6)) will adequately describe viscosity, where t) is viscosity, 

is viscosity at low pressure and fg(p) is a function describing variations in vis- 
cosity. 

8. Pseudo-two- phase thermal properties can be generated by extrapolation of sat- 
uration properties into the near-critical region. 


Heat-Transfer Regions 

1. The near-critical heat-transfer region and its adjacent regions can be delineated 
in an approximate manner on a temperature- entropy state diagram. 

2. It is not possible to set precise near-critical heat-transfer boundaries for the 
following reasons: 

(a) Data are generally insufficient. 

(b) The transition from region to region is not abrupt. 

(c) Most important, the extent of influence of the critical point on heat transfer 
is a function of the process or path by which the fluid approaches the critical point. 

3. Heat transfer in the adjacent regions can be correlated by conventional methods 
appropriate to that region. In fact, the correctable nature of the adjacent regions helps 
determine the near-critical boundaries. 


Near-Critical Heat-Transfer Pecularities 

1. Some characteristics of the near-critical region are 

(a) Nonlinearities in heat flux against temperature difference 

(b) Wall temperature excursions (spikes) 

(c) Similarities to the two-phase regime 

(d) Oscillations 

(e) Large momentum pressure drops 

(f) System- dependent results 

(g) Failure of conventional correlations 

2. Modeling should be required to work throughout region IV, and to exhibit a smooth 
transition to the adjacent regions of figure 8. 
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Free and Natural Convection 


1. Pool results almost universally show an enhancement in heat transfer near the 
critical point. In all cases, the temperature differences are small and the enhancement 
occurs over a limited temperature range and peaks somewhat away from the critical 
point. Large temperature differences diminish the peak. 

2. Free convection systems have no constraining boundaries and consequently re- 
spond favorably to property variations. 

3. Boiling-like behavior in near-critical pools is not accompanied by the marked in- 
crease in the heat-transfer coefficient associated with nucleate boiling. 

4. Most researchers employed conventional techniques away from the critical point. 
However, close to the critical point, modifications such as equations (15) and (16) are 
employed (see table II). 

5. Laminar-to-turbulent oscillations have been observed, but the origin of this phe- 
nomenon has not been explained. These oscillations may well determine the character 
of near-critical oscillations found in natural and forced convection systems. 


Loops - Natural Convection Systems 

1. In the free convection pools and natural convection loops, close proximity of the 
fluid to the critical point appears to enhance heat transfer. This is possibly the result 
of a free boundary system being able to respond to enhancing property variations where 
the temperature differences are low. As will be shown, forced convection systems at 
low temperature differences exhibit similar trends. 

2. Reasonable success has been achieved in correlating the data of both the pools 
and natural convection loops by variable property modifications to the standard equations 
for the given configuration. 

3. The pool experiments have produced considerable support for a pseudo-two- 
phase modeling of a near-critical fluid. On the other hand, the success of single-phase 
correlations indicates either approach may be successful. 

4. To date, no single correlation has been tested against all the available data. 

5. Since these experiments operate very close to the critical point, there is a con- 
tinuing need for improved property data and equations of state. 

6. System oscillations which are commonplace in a natural convection loop near the 
critical point are similar to those observed in forced convection systems. 
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Forced Convection - Experimental 


1. Proximity to the critical point produces a sharp response in wall temperature, 
but the nature of the response is not uniform and the question of increased or decreased 
heat transfer is unclear. 

2. There is evidence of a strong similarity to boiling and there seems to be little 
discontinuity in going from high-pressure boiling to the super- critical region. However, 
since the fluid is not two phase, the formulation is somewhat arbitrary. 

3. Analyses and experiments directed specifically at critical-point heat transfer 
have produced some very interesting variable groupings which deserve continued in- 
vestigation. Among these are 

(a) Integrated specific heat C , and other integrated average properties, 

Jr 


C 


P 


Hw^b 



where H and H, are wall and bulk enthalpy, respectively, and T and T, are wall 
w b wjr.,, w p 

and bulk temperature, respectively. 

(b) Minimum Prandtl number (Pr) m i n 

(c) Eckert reference temperature parameter 


E = 



where T* is the transposed or pseudocritical temperature. 

(d) Turbulence enhancement factors (see the section Theoretical Considerations) 
such as 


dfln p) 


d(ln u + ) 

dy + 

where p is density, y is the coordinate normal to the flow field, y + is dimensionless 
distance, and u + is dimensionless velocity. 

(e) Pseudoboiling parameters such as Stermann number and Martinelli param- 
eter 
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(f) Stability limit (wall temperature spikes) defined by 


q = f(G) (see fig. 20) 

where q is the heat flux and G is mass flow per unit area. 

4. Despite the fact that analysts developing the various correlation or mixing length 
methods have been able to substantiate their work with data, subsequent investigators 
have usually found it necessary to modify the method to fit their own data. No analysis 
is currently widely accepted nor has any analysis been thoroughly tested against the 
available data. 


Geometry Effects 

1. An increase of about 2:1 in the heat-transfer coefficient on the concave side can 
be expected as a result of curvature for gases and fluid systems at high T^/T^. This 
enhancement seems to be a centrifuging effect, particularly in small-bend-angle tubes. 

2. The magnitude of enhancement (or degradation on the convex side) appears to be 
related to the Ito parameter. As a reference for computing enhancement for the dry- 
wall case, correlations of Taylor's type (eq. (9)) appear the most satisfactory. For the 
wet -wall case, the augmentation can shift to the convex wall. 

3. Devices which centrifuge the high-density core to the wall or break up the low- 
density boundary layer, such as twisted tapes and spiral grooves, enhance near-critical 
heat transfer significantly. 

4. Experiment and analysis suggest that a parallel body force should begin to in- 

2 -3 

fluence heat transfer when Gr/Re > 5x10 (where Gr is the Grashoff number and Re 
is the Reynolds number). 

5. Tube entrance effects persist much further down the tube in near- critical fluids 
than in gases. 


Theoretical Considerations in Forced Convection 

1. Conventional mixing length theories must be modified to account for the tremen- 
dous fluid expansion in the near-critical region. These modifications could be in the 
equations, following Goldman's example (ref. 137); or in the eddy diffusivity, following 
Hsu and Smith (ref. 133); or in forming a penetration model, Graham (ref. 150). 

2. Applications of these and other modifications of mixing length theory to near- 
critical data have met with limited success. The theory will probably also have to be 
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made to comprehend upstream history, acceleration, etc. , in order to achieve wide 
success. 

3. Further analysis of accelerated flows should consider the effects of strain rates. 
While the theory of Deissler predicts the proper trends, the magnitudes associated with 
a real near-critical fluid are not correct. Thus, the boundaries and the variation in 
fluid properties must be properly considered (not an easy task). 


Oscillations 

1. Oscillations in near-critical flow can be divided into four groups: laminar to 
turbulent, thermal acoustic, system, and external. 

2. The large density changes of a near-critical fluid are the basic source of the 
oscillating behavior. 

3. The critical frequency for the thermal-acoustic and system oscillations can be 
predicted by utilizing simple concepts of mechanics. 

4. Conventional stability criteria may be employed to predict the onset of thermal- 
acoustic and system oscillations with a fair degree of reliability. 

5. Attempts to predict the amplitude of the oscillations has met with only limited 
success. 


Choking 

1. If a near-critical fluid were in thermodynamic equilibrium, the sonic velocity 
would be low and the flow might be expected to choke very easily. 

2. None of the normal indications of choking have been observed in the near-critical 
experiments. This suggests the possibility of thermodynamic nonequilibrium in these 
flows, which is also consistent with two-phase choking experience. Preliminary nozzle 
data for near-critical nitrogen indicates the expansion to be nonequilibrium. 


Zero-Gravity Operation 

The operation of a 2. 5-liter dewar with internal heaters at supercritical pressure 
and low gravity was reported as not significantly different from operation for earthbound 
tests. Data in this area are lacking. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 15, 1970, 

129-01. 
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APPENDIX A 


SYMBOLS 


A 

constant 

Am 

amplitude parameter, PW n d/Gh^ 

ACT) 

virial coefficient 

A 1 

h 

5' w 

A + 

dimensionless parameter of Van Driest (table V) 

A 

area 

a l 

attraction force constant 

B(T) 

virial coefficient 

b H 

channel half-width 

b l 

repulsion force constant 

C 

constant 

C(T) 

virial coefficient 

C P 

specific heat at constant pressure 

S 

average specific heat, (H w - H b )/(T w - T fe ) 

( c p)„ 

specific heat at constant pressure, 1 atm 

( c p) r 

excess specific heat, C - [C ) 

P \ P /o 

C v 

specific heat at constant volume 

c 

sonic velocity, wave propagation velocity 

D 

self-diffusion coefficient 

D(T) 

virial coefficient 

D ln 

binary diffusion coefficient 

d 

characteristic length (diameter) 

d l 

constant 

E 

Eckert parameter, (T* - T t,)/( T W “ T ^) 


Ep Eg adjustable constants in analytical equation of state (eqs. (4) and (4a)) 
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f 


f() 

fjfa^Cp) 


G 

Gr 


g 

H 

h 

h(Xj) 

h’ 

V** 

h fg 

I 

L 

L’ 

K 


l 

^rp 

m 

Nu 


n 

n l to n 16 
P 


Pr 


P 


* 


energy (heat) eddy diffusivity modification 
momentum eddy diffusitivy modification 
function factor 
functional form 

functions describing varations in thermal conductivity and viscosity, re- 
spectively 

mass flow per unit area, fc>u 

o 9 

Grashoff number: two component, p(p 7 - p)l g/r] ; single component, 

% ATi w 

acceleration 

enthalpy 

heat-transfer coefficient 

function defined by equation (4a) 

modified heat-transfer coefficient, q Q d^‘ G®‘ ® 

latent heat or pseudolatent heat of vaporization 

Ito parameter (eq. (27)) 

axial position, length 

one-half total loop length 

uniform fluid column length between mixing reservoirs 
constant 

characteristic length 
test section length 
mass 

Nusselt number, hd/X 
constant 

constants in modified BWR equation of state 
pressure 

Prandtl number, t^C^/X 

equivalent pressure amplitude when different modes of oscillations occur 
simultaneously 
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Q total heat generated 

q heat flux 

R radius of curvature 

Ra Rayleigh number, GrPr 

Re Reynolds number, pu^d/r] 

(ft gas constant 

r tube radius 

S entropy 

Sr Stermann parameter (boiling number), q Av/u^ AH 
St Stanton number, Nu/RePr 

Sv specific volume number, 

T temperature 

T* transposed critical or pseudocritical temperature (see fig. 6) 
t dimensionless temperature, (T - T c )/T c 

u velocity 

u + dimensionless velocity, u/u* 

u ++ dimensionless velocity (table V) 

u* friction velocity, yv^ 

V volume of gas in resonator cavity 

v specific volume, 1/p 

change in specific volume, fluid to gas 

v* velocity component normal to flow field 

v' average velocity component normal to flow field 
W frequency 

X axial position, length 

X + dimensionless position 

x thermodynamic quality 

x Q adjustable constant in analytical equation of state (eqs. (4) and (4a)) 
x 


dimensionless parameter of equation (4), 


P/Pj 
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y 


coordinate normal to flow field 


y + 

y ++ 

z 

Z 

a 

p 

h 

r 

+ 

y 

A() 

5 



i H(P, T) 

Mi 


dimensionless distance, u*y/u Q 
dimensionless distance (table V) 
compressibility factor, P/pdTT 
pressure head between heater and cooler 
void fraction (two-phase system) 
reference temperature, Q 0 u */^pj t o T o 


expansion coefficient 

adjustable constant in analytical equation of state (eqs . (4) and (4a)) 

ratio of specific heats (isentropic expansion coefficient), Cp/C v 

changeover point in equation (42) 

change in a function 

boundary- layer thickness 

adjustable constant 

dimensionless boundary- layer thickness, u + 6/v 

momentum eddy diffusivity 

energy (heat) eddy diffusivity 

dimensionless eddy diffusivity, e/v 

dimensionless eddy diffusivity (Deissler form, eq. (42)) 

dimensionless eddy diffusivity (Van Driest form, eq. (41)) 

viscosity 

viscosity (p — 0, T) 
thermal conductivity 

thermal conductivity given by equation (5) 
dissociation contribution to thermal conductivity 
thermal conductivity (p — 0, T) 
chemical potential 

Joule-Thompson coefficient, = (dT/9P)jj = 0 
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V 

4 

p 

T 

T contact 

<P 


kinematic viscosity, rj/p 
dimensionless distance, §A>jj 
density 
shear force 

fluid residence time (eq. (45b)) 
property parameter (eq. (24)) 
integral averaged property 


% 

(XfP' 

d) 

Subscripts: 

av 

b 

c 

calc 

core 

4 

cr 

cv 

cx 

f 

fm 

H 

g 

l 

min 

n 

o 


Martinelli parameter, turbulent-turbulent flow; 


1 - x 


0 - 9 /P2\°->i\ 0 - 1 


\P b 




modified Martinelli parameter, where x is pseudoquality 
mass flow rate 


average 

properties evaluated at bulk conditions 

critical 

calculated 

core of fluid 

centerline 

curved 

concave surface 
convex surface 

properties evaluated at film conditions 
properties evaluated at mean film conditions 
at constant enthalpy 
gas 

properties evaluated at liquid conditions 
properties evaluated according to equation (22a) 
natural 

properties evaluated at wall conditions 
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p penetration 

pg perfect gas 

sat saturation 

st straight 

test conditions of experimental test 

tot total 

w properties evaluated at wall conditions 

x properties evaluated at reference conditions 

5 edge of boundary layer 

cp properties evaluated at void weighted conditions 

0 initial 

1 heavy species I Do not apply to equations (4), (4a), (5a), and (6a), 

2 light species J where 1 and 2 are ordinary subscripts. 

°° free stream 
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APPENDIX B 


HEAT-TRANSFER REGIMES AND THE PROCESS MAP 

The heat-transfer problems can be loosely classified into the four regimes, as 
illustrated in figure 8. However, it is not sufficient just to know what regime is of in- 
terest, the process path within these regimes must also be known. For example, a 
process in region IV (fig. 8) may exhibit very few of the pecularities discussed in the 
text if the process traverses only a small part of this regime. As a consequence, the 
designer would have little difficulty in applying a correction to the ordinary Dittus- 
Boelter equation (ref. 138) and in predicting the heat transfer with reasonable accuracy. 

Consider first the problem of establishing the regime of interest, and second the 
problem of establishing the process path. 

In establishing the regime of interest, usually only the thermodynamic limits must 
be examined (i. e. , the operating temperature, the pressure, and the enthalpy (or en- 
tropy rise). With these data and a figure similar to figure 8, some decision can be 
made as to the regime. Such a procedure is flow charted as figure 35. Basically, the 
message of figure 35 is that if you can avoid operating in the near-critical regime you 
are better off as far as predicting heat transfer is concerned; however, if you must op- 
erate in this regime, the material of this text should be beneficial in delineating prob- 
lems which must be considered. 

In determining the process path assume that a figure similar to figure 8 has been 
constructed, for example, figure 36. Now examine the five typical process paths illu- 
strated in the figure: 

(1) Process path 1: This represents the standard case of a gas entering the appara- 
tus (perhaps with moderate precooling) and exiting as a low-Mach-number gas. The 
heat-transfer and pressure-drop equations are available and cited herein (see the main 
text, the section HEAT-TRANSFER REGIONS). 

(2) Process path 2: Such a process is closely related to the standard case, except 
that the working fluid is a liquid at the system inlet and a subcooled liquid at the system 
exit. This fluid may be at either subcritical or supercritical pressure in a system with 
low temperature differences - and, of course, there would be no change in phase or 
pseudophase. Again, standard equations similar to those presented in the main text 
(see sections HEAT-TRANSFER REGIONS and the NEAR CRITICAL HEAT-TRANSFER 
REGION) should apply to yield heat transfer and pressure drop. 

(3) Process path 3: In such a process, there is a change of phase which complicates 
the prediction of heat transfer and pressure drop. The process starts as a subcritical, 
subcooled liquid and terminates in the two- phase regime. While it is difficult to assess 
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Figure 35. - Flow chart showing procedure of characterizing the heat -transfer process within the four heat-transfer regimes. 




















Process 1*. Starts and ends as a gas away from critical (sub- or supercritical pressure) 

Process 2: Starts and ends as a liquid away from critical (sub- or supercritical pressure - low AT) 

Process 3, 3a : Starts liquid or two phase and ends two phase or gas 

Process 4: T i n |et <T * <T outlet and hi 9 h AT and hi 9 h heat fluxes. Usually starts as near-critical liquid 




(b) Superposition of process on region map. 


Figure 36. - Typical process paths and characteristic heat-transfer regions. 


the region around the saturation locus, several analyses, as discussed in the text, are 
available which predict heat transfer and pressure drop with reasonable accuracy. 

(4) Process path 3a: Such a process enters the system as a two- phase fluid; how- 
ever, the exhausted fluid is gaseous or perhaps droplet- laden. This is a difficult regime 
to assess; however, the techniques of Hsu, et al. (ref. 49) should be applicable. 

(5) Process path 4: These are the systems operating in the near-critical regime, 
where, as a general rule, heat-transfer predictions are poor. However, if the bound- 
aries of region IV are not crossed (see fig. 36), it appears to be possible to make a 
reasonable prediction. If in the final analysis, the process crosses the boundaries, the 
problem will be a real challenge. Some of the problems involved in analyzing such a 
system are covered in the main text of this report. 
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